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Preface 
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helped  type  drafts  of  the  thesis  but  also  endured  six  months 
of  a grumpy,  nervous,  and  seldom  seen  husband. 
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Abstract 

An  investigation  of  the  forecasts  of  aeronautical  price 
indices  done  by  the  Aeronautical  Systems  Division  (ASD)  of 
the  United  States  Air  Force  Systems  Command  was  performed. 

The  state  of  the  art  in  correction  for  autocorrelation  was 
summarized.  The  historical  price  indices  used  by  ASD  were 
recalculated.  Using  the  recalculated  historical  composite 
price  indices,  different  forecasting  techniques  were  examined 
including:  simple  linear  regression,  multivariate  regression, 

regression  with  nonlinear  variables,  Generalized  Least 
Squares  with  the  first  order  autoregressive  model,  first  dif- 
ferences, and  per  cent  change  regression.  The  Gross  National 
Product  price  deflator,  the  Primary  Metals  Index,  the  Elec- 
trical Machinery  Index,  the  national  total  unemployment  rate, 
and  the  money  supply , Ml , were  all  tried  as  independent  var- 
iables in  the  regressions.  The  Durbin  Two  Step,  the  Search 
technique,  and  a Bayesian  technique  were  all  used  to  estimate 
the  first  order  autoregressive  coefficient,  p.  When  start- 
ing with  the  composite  price  indices,  the  per  cent  change  re- 
gression, with  the  Gross  National  Product  price  deflator  as 
the  only  independent  variable,  was  recommended  as  the  best 
forecasting  technique.  Prediction  intervals  were  calculated. 
The  effect  of  errors  in  the  Wharton  Econometric  Forecasting 
Associates  prediction  of  the  Gross  National  Product  price  de- 
flator was  discussed. 
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AN  INVESTIGATION  OF  THE  FORECASTING  OF 
AERONAUTICAL  PRICE  INDICES 

I . Introduction 

In  1973  the  Aeronautical  Systems  Division  (ASD)  of  the 
United  States  Air  Force  Systems  Command  begin  publishing  an- 
nual reports  of  forecasted  aeronautical  price  indices.  In- 
cluded in  the  reports  were  airframe,  engine  and  avionics  in- 
dices for  both  production  and  development.  These  forecasted 
indices  were  used  by  budget  estimators  in  ASD  to  convert  bud- 
get estimates  in  current  dollars  to  budget  estimates  in  then- 
year  dollars.  For  example,  a budget  estimator  in  an  ASD  Sys- 
tem Program  Office  might  estimate  that  the  budget  for  a spe- 
cific aircraft  avionics  production  program  would  be  $100,000 
in  1976  dollars.  If  one  assumes  that  a forecasted  avionics 
production  price  index  for  1977  was  110  with  a base  of  100 
in  year  1976,  then  a budget  estimator  would  submit  a budget 
for  avionics  production  of  $110,000  in  1977  dollars.  This 
budget  would  then  be  sent  to  higher  levels  for  approval. 

This  thesis  deals  directly  with  the  ASD  reports  of  fore- 
casted price  indices.  The  thesis  covers  two  major  areas: 
one,  the  accuracy  of  the  predictions,  and  two,  the  possibil- 
ity of  improved  accuracy  through  different  or  advanced  sta- 
tistical techniques.  The  accuracy  is  checked  by  starting  at 
the  beginning  with  a recalculation  of  the  historical  indices, 
and  then  working  all  the  way  through  to  an  estimate  of  the 
prediction  intervals . Different  statistical  techniques  are 
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examined  including  the  addition  of  extra  explanatory  vari- 
ables and  the  assumption  of  a first  order  autoregressive 
model.  All  mathematical  techniques  are  clearly  explained 
and  important  computer  programs  are  listed. 
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I I . Background 

Aeronautical  Systesm  Division  (ASD)  has  published  annual 
reports  of  forecasted  price  indices  since  1973.  The  first 
two  reports,  denoted  110A  and  HOB,  used  a similar  methodol- 
ogy. The  HOC  report  used  a slightly  different  technique 
(Ref  1;  Ref  2;  Ref  3) . 

110A  and  110B  Reports 

ASD  used  a straightforward  three  step  methodology  in 
their  first  two  reports,  110A  and  110B.  First,  they  gathered 
historical  data  and  developed  six  price  indices.  Second, 
they  regressed  these  six  price  indices  against  the  Gross 
National  Product  (GNP)  price  deflator.  Third,  they  used  a 
prediction  of  the  GNP  price  deflator  to  forecast  the  price 
indices 

In  step  one,  ASD  used  historical  data  from  aerospace 
contractors,  combined  with  nationally  published  statistics. 

To  develop  the  six  composite  historical  indices,  Production 
Engine,  Production  Airframe,  Production  Avionics,  Develop- 
ment Engine,  Development  Airframe,  and  Development  Avionics, 
ASD  used  a weighted  average  of  subindices.  For  example,  to 
construct  the  Development  Engine  index  ASD  used  the  follow- 
ing procedure.  First ; they  constructed  a purchased  parts 
subindex  using  the  following  formula: 

(.6)  Engine  + (.4)  Raw  Materials  = Purchased  Parts  (1) 
Labor  Subindex  Subindex 

Subindex 
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Second,  they  used  this  purchased  parts  subindex  along  with 
other  subindices  to  form  the  Development  Engine  index. 

(.4484)  Engine  + (.1632)  Engineering 

Labor  Subindex  Labor  Subindex 

+ (.1942)  Raw  Materials  + (.1942)  Purchased 
Subindex  Parts 

= Development 

Engine  Index  (2) 

The  subindices  were  developed  using  data  from  the  Bureau  of 
Labor  Statistics  and  industry.  Engine  labor  data  were  from 
Bureau  of  Labor  Statistics  code  3722.  The  Raw  materials  in- 
dex was  derived  from  several  wholesale  price  index  categories. 
Engineering  labor  data  was  obtained  from  aerospace  contrac- 
tor data.  The  weights  and  methodology  used  were  based  on 
contractor  data,  and  on  a Rand  Corporation  report  published 
in  1970  (Ref  4) . 

In  step  two,  ASD  related  the  historical  data  to  the 
historical  Gross  National  Product  price  deflator  (GNP^) . 

They  did  this  by  using  simple  linear  regression.  ASD  then 
had  six  regression  equations,  engine,  airframe,  and  avionics 
for  both  production  and  development.  The  results  for  the  De- 
velopment Engine  index  were: 

Development  Engine  Index  = -10.8  + 1.073  GNP^  (3) 

R , the  coefficient  of  determination,  was  .982.  Using  the 
criterion  of  R , it  appeared  that  the  regression  had  resulted 
in  a very  good  fit  of  the  data.  Actually,  as  will  be  seen 
later,  the  data  was  autocorrelated  and  the  fit  was  not  nearly 
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2 

as  good  as  was  first  indicated  by  the  high  R . 

In  step  three,  ASD  attempted  to  find  a reliable  predic- 
tor of  the  GNP  deflator  which  could  be  inserted  into  the  re- 
gression equations  in  order  to  predict  the  appropriate  index. 
ASD  chose  the  Wharton  Econometric  Forecasting  Associates, 

Inc.  prediction  of  the  GNP  deflator.  ASD  did  not  use  Whar- 
ton's raw  predictor,  but  smoothed  the  predictor  before  using 
it.  ASD  smoothed  the  predictor  by  fitting  an  exponential 
curve  to  it.  ASD  said  they  did  this  to  reduce  the  impact  of 
short  term  fluctuations  (Ref  2:47).  That  then  is  the  basic 
technique  used  by  ASD  in  their  first  two  reports,  110A  and 
HOB. 

110C  Report 

The  110C  report,  published  in  1976,  had  three  major 
changes.  The  first  change  was  that  ASD  regressed  where  pos- 
sible each  subindex  against  what  appeared  to  be  an  appropriate 
independent  variable.  Three  independent  variables  were  used: 
GNP  price  deflator,  Primary  Metals  Price  Index,  and  the  Elec- 
trical Machinery  Price  Index.  The  second  major  change  was 
that  a per  cent  change  regression  was  used  in  the  attempt  to 
overcome  apparent  autocorrelation.  Per  cent  changes  in  each 
dependent  variable  were  regressed  against  per  cent  changes, 
in  the  independent  variable.  Third  and  last,  the  unsmoothed 
Wharton  GNP  deflator  was  used.  As  a result,  instead  of 
three  major  indices  regressed  for  both  production  and  develop- 
ment for  a total  of  six  regression  equations,  all  but  the 
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overhead  subindices  were  regressed  for  a total  of  nine  regres- 
sion equations.  These  regression  equations  were  used  to  pre- 
dict future  values  of  the  subindices.  The  overhead  subin- 
dices were  predicted  using  the  historic  compound  rate.  For 
example,  the  historical  7.3%  Overhead  subindex  was  calculated 
by  compounding  7.3%  back  and  forth  from  a 1967  base  value  of 
100.  The  7.3%  Overhead  subindex  for  1975  was 

(100)  (1  + . 073) 8 = 175.71  (4) 

In  the  HOC  report,  this  same  compound  rate  was  used  to  pre- 
dict future  values  of  the  index.  Thus  the  7.3%  subindex  for 
1976  was 

(100)  (1  + . 073) 9 = 188.54  (5) 


Weights  were  applied  to  the  subindex  predictions  to  develop 
the  six  composite  indices.  The  methodology  of  the  110A, 
110B  and  110C  reports  is  contrasted  in  Fig.  1. 


Future  Interests 


The  per  cent  change  regression  used  in  report  110C  ap- 
parently corrected  the  data  for  autocorrelation.  ASD,  how- 
ever, has  expressed  an  interest  in  other  regression  techni- 
ques that  are  simple,  that  correct  for  autocorrelation,  and 
that  are  more  effective  than  the  per  cent  change  regression. 

In  the  110C  report,  the  subindex  was  first  regressed, 
the  predictions  were  made,  and  then  the  predictions  were 
weighted  and  aggregated  to  form  the  composite  indices.  In 
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Methodology 


110A  and  HOB  110C 


1.  Form  historical  subin-  1.  Form  historical  subin- 

dices.' Weight  these  subin-  dices. 

dices  to  form  six  composite 
indices . 

2.  Regress  each  of  the  six  2.  Using  per  cent  changes, 

composite  indices  against  the  regress  all  of  the  subindices, 
historical  GNP  price  defla-  except  overhead  and  material, 
to”.  against  the  GNP  price  deflator. 

Regress  the  materials  subin- 
dices against  either  the  Pri- 
mary Metals  index,  or  the  Elec- 
trical Machinery  index. 

3.  Use  the  smoothed  predic-  3.  Project  the  overhead  sub- 
tions  of  the  GNP  price  defla-  indices  into  the  future  using 
tor  in  the  regression  equa-  the  historical  compound  rate, 
tions  to  predice  the  future 

values  of  all  six  indices. 

4.  Use  unsmoothed  predictions 
of  the  independent  variables 
in  the  regression  equations 

to  predict  per  cent  changes  in 
the  subindices. 

5.  Change  per  cent  changes 
into  subindices. 


6.  Weight  subindices  to  pro- 
duce the  predicted  values  for 
the  six  composite  indices. 
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the  previous  reports,  the  subindices  were  first  aggregated 
and  then  the  resulting  composite  indices  were  regressed. 

ASD  has  expressed  an  interest  in  the  comparison  of  two  tech- 
niques, both  using  per  cent  change  regressions.  The  first 
technique  would  be  to  predict  then  aggregate,  and  the  second 
technique  would  be  to  aggregate  then  predict. 

This  thesis  will  investigate  two  areas  in  which  ASD  has 
expressed  interest: 

1)  Are  there  other  techniques  that  can  and  should  be 
used  to  correct  for  autocorrelation? 

2)  Which  is  better,  regress,  then  aggregate,  or  aggre- 
gate, and  then  regress? 

Thesis  Outline 

A recalculation  of  the  historical  indices  was  done  in 
Chapter  IV.  In  Chapter  V the  per  cent  change  regression  as 
done  in  the  HOC  report  was  repeated  with  one  addition,  pre- 
diction intervals  were  given.  Also  in  Chapter  V,  two  dif- 
ferent approaches  to  per  cent  change  regression  were  dis- 
cussed. The  first  approach  was  to  regress,  and  then  aggre- 
gate the  predictions.  The  second  approach  was  to  aggregate 
the  data,  and  then  regress.  Chapter  VI  investigated  other 
regression  techniques.  Chapter  VII  discussed  the  error  in 
the  Wharton  predictors  and  its  possible  effect  on  the  error 
of  the  predicted  indices.  Finally  Chapter  VIII  summarized, 
presented  significant  findings,  and  made  recommendations  for 
future  study. 
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III.  Regression  Theory 


Preface 


This  chapter  has  been  prepared  in  order  to  furnish  a 
concise,  useful,  and  accurate  summary  of  applicable  mathe- 
matical theory.  This  chapter  first  provides  a discussion  of 
the  assumptions  made  when  doing  "ordinary",  regression.  It 
then  provides  a detailed  discussion  of  autocorrelation,  a 
violation  of  one  of  the  assumptions  of  regression  analysis. 


Assumptions  of  the  General  Linear  Model 


The  assumptions  of  the  General  Linear  Model  (sometimes 
known  as  Ordinary  Least  Squares,  OLS)  can  be  placed  in  two 
sets.  The  first  set  allows  for  efficient  estimates.  The 
second  set,  which  contains  only  one  assumption,  allows  for 
easy  computation  of  statistical  tests  on  the  model  and  on 
its  predictions.  The  first  set  of  assumptions,  along  with 
brief  explanations,  are  listed  below: 

^nxl  - ^nxk^kxl  + £nxl 

Y is  expressed  as  a linear  combination  of  k,  B's  with  an  er- 
ror term,  ie.,  the  model  is  correctly  specified. 


2>  E<£nxllW  ■ °nxl 


The  expected  value  of  the  error  is  zero. 


3>  V<EnxllXnxk>  • 01 


* X 


r 


The  errors  are  independent  of  each  other  and  all  have  the 
same  variance. 


4)  has  rank  k 

The  columns  of  X are  independent,  ie.,  the  X variables  are 


independent  of  each  other. 


5)  may  consist  of  either  fixed  or  random  variables 

where  A^  stands  for  a matrix  A,  b rows  by  c columns  (Ref  5: 

111) . 

When  these  assumptions  are  satisfied.  Ordinary  Least 
Squares  will  yield  the  best  linear  unbiased  estimates  of  the 
coefficients,  B^xi  (Ref  5:119).  One  other  assumption  is  fre- 
quently made. 

6)  Probability  CY ! Xnxic®icxi ) normally  distributed. 

The  errors  are  normally  distributed. 

This  assumption  means  that  when  OLS  is  performed,  estimates 
of  the  coefficients,  B^xi*  will  be  maximum  likelihood  esti- 
mates. Also  confidence  intervals,  prediction  intervals,  and 
statistical  tests  can  be  calculated. 

Checks  on  the  Assumptions  of  Ordinary  Least  Squares 

In  this  thesis,  both  visual  and  statistical  tests  were 
performed  to  check  on  the  applicability  of  the  six  assump- 
tions or  Ordinary  Least  Squares.  OMNITAB,  the  computer  pro- 
gram language  used  in  this  thesis,  provided  four  plots  that 
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were  used  to  visually  check  assumptions  1,  3 and  6 (Ref  6: 
138-146).  The  four  plots  were  1)  residuals  versus  time,  2) 
residuals  versus  the  predicted  values,  3)  residuals  versus 
the  independent  variable,  and  4)  a probability  plot  of  the 
residuals.  These  plots  were  utilized  in  the  following  way. 

If  the  points  on  plot  1 were  not  scattered  around  zero,  then 
the  errors  may  not  have  been  independent  and  assumption  3 
was  assumed  to  be  violated.  If  the  points  on  plot  2 were 
not  scattered  around  zero,  assumption  1 was  assumed  to  be 
violated,  an  indication  that  a different  model  such  as  a 
curvilinear  regression  might  have  done  a better  job.  If  the 
points  on  plot  3 were  not  scattered  around  zero,  the  vari- 
ances of  the  errors  may  not  have  been  constant,  and  assump- 
tion 3 was  assumed  to  be  violated.  If  the  points  on  plot  4 
did  not  lie  on  a straight  line,  then  assumption  6 was  assumed 
to  be  violated,  an  indication  that  the  errors  may  not  have 
been  normally  distributed. 

In  addition  to  these  visual  tests,  numerical  tests  were 
also  used  to  check  on  the  assumptions.  A t test  was  used  to 
test  the  significance  of  the  coefficients  in  the  regression 
and  to  see  if  assumption  1 may  have  been  violated  (Ref  5:138). 
The  OMNITAB  program  also  provided  estimates  of  the  accuracy 
of  its  calculations,  which  were  useful  when  assessing  the  in- 
dependence of  the  columns  of  the  X matrix,  assumption  4 (Ref 
6:149) . 


The  last  numerical  test  used  the  Durbin  Watson  statis- 
tic (Ref  5:199-201).  The  test  checks  to  see  if  the  errors 
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are  not  independent,  but  rather  are  related  to  each  other 
in  time,  a violation  of  assumption  3.  If  the  error  terms 
are  positively  related  to  each  other  in  time,  then  they  are 
said  to  be  positively  autocorrelated . If  the  error  terms 
are  negatively  related  to  each  other  in  time,  they  are  said 
to  be  negatively  autocorrelated.  This  last  numerical  test 
is  a statistical  test  consisting  of  a null  hypothesis,  Hq, 
an  alternative  hypothesis,  H^,  a test  statistic,  D,  and  a 
rule.  The  test  for  positive  autocorrelation  is: 


Hn  : No  positive  autocorrelation 


: Positive  autocorrelation 


D = 


n-1 

£.  <ea*l 
a=l 


- e ) 
cr 


n 


a=l 


a 


Rule  if  D < 

if  DL  < D < ^ 


if  D > D 


U 


Reject  Hq 
Inconclusive 
Do  not  reject  HQ 


(6) 


where  e^  are  the  residuals  from  Ordinary  Least  Squares.  The 
test  for  negative  autocorrelation  is  the  same  with  4-D  re- 
placing D.  Dl  and  were  computed  and  tabled  by  Durbin  and 
Watson. 

It  is  important  to  realize  that  when  D > D^  some  authors 
state  that  the  user  should  accept  Hq  (Ref  7:358).  Strictly 
speaking,  this  is  not  correct.  There  is  little  one  can  say 
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statistically  when  D > D^.  The  significance  levels  given  in 
the  D^,Djj  table  apply  only  to  the  rejection  of  Hq.  In  prac- 
tice, however,  when  D approaches  2.0  the  user  can  be  fairly 
confident  that  the  errors  are  not  autocorrelated . 

Assumption  2 was  indirectly  tested  when  the  assumption 
of  the  correctness  of  the  model  was  tested.  Assumption  5 
was  not  tested.  The  historical  independent  variables  were 
assumed  fixed.  The  treatment  of  the  error  in  the  predicted 
values  of  the  independent  variables  is  included  in  Chapter 
VII. 


What  to  do  When  Autocorrelation  is  Detected 

There  are  many  things  the  analyst  can  do  when  autocor- 
relation is  detected.  He  can  ignore  it  and  perform  Ordinary 
Least  Squares  as  usual,  or  he  can  check  to  see  if  the  origi- 
nal model  is  misspecif ied , or  he  can  assume  a first  order 
autoregressive  model.  If  the  user  assumes  a first  order 
autoregressive  model,  there  are  many  techniques  available 
for  use.  All  of  these  approaches  are  discussed  below. 

Ignore 

The  analyst  can  ignore  the  warning  of  possible  autocor- 
relation and  perform  Ordinary  Least  Squares  on  the  data. 

This  will  yield  estimates  of  the  coefficients;  but,  if  the 
data  is  autocorrelated,  the  estimates  may  not  be  the  best 
estimates  available.  Also  if  the  data  is  autocorrelated, 
Ordinary  Least  Squares  will  yield  calculated  variances  that 
underestimate  the  true  variances  (Ref  5:254-257).  Fig.  2 
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is  an  example  of  what  might  happen  if  Ordinary  Least  Squares 
is  performed  on  autocorrelated  data. 

Misspecif ied  Model 

A low  Durbin  Watson  test  statistic  indicating  autocor- 
relation may  result  from  a misspecif ication  of  the  model.  A 
variable  may  be  missing  or  a variable  may  be  in  the  wrong 

form.  For  example  if  the  real  relationship  that  exists  be- 

2 

tween  Y and  X is  Y = X , and  a straight  line  is  fitted  to 
the  data,  and  if  time  is  highly  related  to  X,  then  the  data 
may  appear  to  be  positively  autocorrelated.  The  real  prob- 
lem is  not  autocorrelation,  but  rather  a misspecification 
of  the  model.  See  Fig.  3. 

Assume • First  Order  Autoregressive  Model 

If  one  is  not  satisfied  with  ignoring  the  detection  of 
possible  autocorrelation,  and  if  different  models  and  addi- 
tional variables  have  been  tested  to  no  avail,  then  one  recom- 
mended next  step  is  to  assume  a first  order  autoregressive 
model.  The  first  order  autoregressive  model  is: 

^nxl  ~ ^nxk^kxl  + £nxl 
£t  = p£t-l  + yt 

where  y^  is  a normally  distributed  independent  random  vari- 
able with  zero  mean  and  constant  variance,  e is  the  error 
at  time  t (Ref  5:250).  In  other  words,  the  error  this  peri- 
od is  equal  to  a constant  times  the  error  last  period  plus  a 
random  error.  Usually  another  assumption  is  included  in  the 


(7) 

(8) 
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first  order  autoregressive  model.  It  is  that  | p | < i (Ref  5: 
250).  In  this  thesis,  however,  no  restriction  will  be  put 
on  p.  Techniques  will  be  discussed  that  assume  p = 1, 

| p | < 1 , and  -oo<p<». 


| p | < 1 Assumption 

The  assumption  of  | p | < 1 allows  for  cases  of  first  order 

autocorrelation  where  the  error  does  not  go  to  infinity  in 

time.  Without  going  into  a lot  of  detail,  it  is  sufficient 

to  say  that  when  p is  known  a matrix  called  P can  be  con- 
' nxn 

structed  (Ref  5:252-254). 


P 

nxn 


/ITfF  0 0 

-p  1 0 

0 -Pi 


0 0 
0 0 
0 0 


0 


0 


0 


-p 


1 


where  |p|<l.  This  matrix  can  be  used  to  transform  the  orig- 
inal equation,  Eq  (7),  into  a new  equation  with  independent 
error  terms. 


P Y . 
nxn  nxl 


P X , + 

nxn  nxk  kxl 


P e 
nxn  nxl 


(9) 


or  for  a model  with  a constant  and  a linear  term 
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“ — 

m mm 

p — * 

/r^rY1 

y2-pyi 

/I -Q* 
1-P 

60 

/1-Pix1 

y2-pXi 

*1 

/r^p^iy 

e2‘pel 

VpYn-l 

1 

1-p 

2 

+ 

VpXn-l 

3 

+ 

_en‘pen-i 

(10) 


Applying  Ordinary  Least  Squares  to  the  transformed  data,  the 
data  in  the  braces  labeled  1,  2,  3,  yields  the  best  estimate 
of  and  8-^.  This  procedure  is  a special  case  of  General- 
ized Least  Squares.  The  problem,  of  course,  is  that  p is 
not  normally  known.  Various  techniques  have  been  invented 
to  estimate  both  p and  - Four  of  these  techniques  will 

be  discussed  here:  Cochrane-Orcutt , Durbin  Nonlinear,  Durbin 

Two  Step  and  a method  called  Search. 

Cochrane-Orcutt . One  technique  discussed  by  Cochrane 
and  Orcutt,  but  not  recommended  by  them  because  of  possible 
bias,  was  to  start  with  a guess  of  p in  Eq  (9)  and  then  esti- 
mate by  Ordinary  Least  Squares  (OLS)  (Ref  9:53).  Then 

the  estimate  of  was  put  back  in  the  original  equation, 

Eq  (7),  and  the  resulting  residuals  were  used  to  estimate  p. 
This  new  p was  then  put  back  in  Eq  (9)  and  a new  was 

estimated  using  OLS.  The  process  was  continued  until  suc- 
cessive estimates  of  p differed  by  a small  amount.  This 
method  has  been  frequently  illustrated  without  the  first  row 
of  Eq  (9).  If  the  autocorrelation  process  has  been  going  on 
before  the  first  observation  point,  then  leaving  out  the 
first  row  reduced  the  accuracy  of  this  technique  (Ref  10:96). 
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Durbin  Nonlinear . Two  of  the  other  techniques  discussed 
here  that  estimated  p and  8^^  also  do  not  use  row  1 of  Eq 
(10).  One  technique  is  called  the  Durbin  Nonlinear  (Ref  11). 
For  the  case  of  a constant  plus  a linear  term,  Durbin  rewrites 
Eq  (9)  as 

Yt  - e0Cl-P)*B1Xt  . »Yt.1  . ut  (11) 


Yt  - Bq ' . BlXt  • B2Xt.1  . BjY^j  * «t 


where  y is  a random  error.  The  Nonlinear  technique  mini- 
mizes the  sum  of  squares  of  the  errors  in  Eq  (12)  with  the 
restriction  that  8 -,  / 6 “83.  63  will  be  the  estimate  of  p. 

If  desired,  this  p can  be  used  in  Eq  (9)  to  estimate  8^x^. 

This  technique,  essentially  least  squares  subject  to  a quad- 
ratic constraint,  was  tried  with  little  success  by  Potluri 
and  Griliches  (Ref  11).  This  technique  is  difficult  to  pro- 
gram on  a computer. 

Durbin  Two  Step . The  other  technique  discussed  here 
that  estimates  p and  8^x^  without  row  1 is  called  the  Durbin 
Two  step  (Ref  12:150-153).  Durbin  recommends  directly  per- 
forming least  squares  on  Eq  (12)  while  ignoring  any  constraint 
on  the  coefficients.  The  estimate  of  6^  is  then  used  as  the 
estimate  of  p.  If  desired  this  p can  be  used  in  Eq  (9)  to 
estimate  6^x^.  The  Durbin  Two  step  is  certainly  simple  and 
is  more  effective  than  it  first  may  appear.  It  leads  to  es- 
timates having  the  same  efficiency  as  the  least  square  esti- 
mates in  large  samples  (Ref  11;  Ref  12:150-153). 
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Search . The  last  technique  discussed  under  the  assump- 
tion that  | p | < 1 is  called  the  Search  technique.  Like  the 
Durbin  Nonlinear  and  the  Durbin  Two  step  it  attempts  to  mini- 
mize the  sum  of  squares  of  error  of  Eq  (9) , but  unlike  the 
Durbin  techniques  it  includes  row  one  of  Eq  (9) . The  Search 
technique  is  accurate.  The  minimum  sum  of  squares  is  always 
found.  Most  important  of  all,  the  Search  technique  yields 
the  maximum  likelihood  estimates  of  p and  6^*1  (Ref  13:  Ref 
14).  The  Search  technique  is  simple  to  understand.  For  suc- 
cessive increments  of  p between  1 and  -1,  excluding  0,  OLS 

is  performed  on  Eq  (9)  yielding  estimates  of  p,  and 

2 

and  the  variance,  S . There  will  be  only  one  minimum  point 
2 

for  S . p,  and  corresponding  to  this  minimum  will  be  the 

maximum  likelihood  estimates  (Ref  13) . 

The  Search  technique  will  probably  require  more  itera- 
tions than  Cochrane-Orcutt  or  the  Durbin  techniques;  however, 
with  careful  planning  the  number  of  iterations  can  be  re- 
duced. If  one  knows  the  sign  of P (ie.,  positive  or  negative 

autocorrelation)  then  only  one  half  the  interval  of  p need 

2 

be  searched.  Since < there  is  only  one  minimum  S point,  not 

every  increment  of  p need  to  evaluated.  The  interval  of 

under  search  can  be  bisected  and  the  midpoint  and  endpoints 

2 

can  be  evaluated  for  S . The  two  resulting  intervals  can  be 

2 

bisected  and  their  midpoints  can  be  evaluated  for  S . The 

2 

process  continues  and  the  minimum  point  of  S is  located  in 
successively  smaller  intervals.  Knowing  the  sign  of  p and 
using  the  above  process  can  reduce  the  maximum  number  of  it- 
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erations  needed  to  estimate  p.  No  more  than  14  iterations 

are  necessary  to  estimate  p to  2 decimal  places. 

Prediction.  Goldberger  showed  how  to  predict  future 

2 

values  of  Y along  with  their  prediction  variances,  using 
p and  determined  from  the  previous  four  techniques  (Ref 

15:372-374).  The  best  linear  unbiased  estimate  of  future 
values  of  Y is  given  by 


Yt  ' Xlxk,t  Bkxl  * <Y„  - XnBkxl> 


(13) 


th 


where  X^xk  t are  the  predictor  variables  for  the  t period 


in  the  future,  and  X is  the  last  row  of  X 


n 


nxk ' 


Y is  the 
n 


th 


last  row  of  Y^^,  and  Y^  is  the  predicted  value  for  the  t 

2 

period.  The  variance  of  this  prediction,  at  , is  given  by 


<Jt2  = (1  - p2t)a2  + (1  - p2)a2  [1  + (Xlxk>t  - pXnK 

(14) 

(Xnxk  pnxn  Pnxn  Xnxk)  1(Xlxk,t  ' pXn5 1 

2 2 

where  Pnxn  is  as  given  in  Eq  (9),  and  (1  - p )a  is  the  vari- 
ance of  the  errors  in  Eq  (9) . 

This  concludes  the  discussion  of  the  various  techniques 
used  if  one  decides  to  use  the  first  order  autoregressive 
model  with  the  assumption  that  | p | < 1 . Next  a technique  used 
to  estimate  p for  the  assumption  that  -®<p<«>  will  be  discus- 
sed. 


— oo< p < oo  Assumption 

The  assumption  that  -<»<p<oo  is  rarely  discussed  in  the 
literature,  p greater  than  one  means  the  errors  go  to  in- 
finity in  time.  This  is  sometimes  called  an  explosive  pro- 
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cess.  Although  an  explosive  process  may  seem  to  be  unlikely, 
some  models  may  have  errors  that  increase  slightly  with  time. 

Bayesian.  Bayesian  techniques  have  been  used  to  develop 
the  probability  of  p given  the  data  for  the  assumption  that 
-oo<p«x>  (Ref  16)  . The  result  is  that  the  probability  of  p 
given  the  data,  is  inversely  related  to  the  variance  of  the 
error  resulting  from  OLS  performed  on  Eq  (9)  without  row  one. 
The  smallest  variance  will  correspond  to  the  most  likely  es- 
timate for  p.  Unfortunately,  if  this  estimate  of  p is  great- 
er than  +1  or  less  than  -1,  then  the  theory  stops.  No  one 
has  discussed  how  to  predict  future  values  if  | p | > 1 . How- 
ever, the  Bayesian  technique  was  used  in  this  thesis  to  get 
an  idea  of  how  big  p was  when  other  estimates  of  p indicated 
p was  greater  than  .99. 

p = 1 Assumption 

The  simplest  of  all  techniques  used  to  correct  for  auto- 
correlation is  also  the  most  restrictive.  It  assumes  p = 1. 
The  technique  is  called  first  differences.  For  illustration 
assume  a model  with  a constant  and  a linear  term.  Subtract- 
ing successive  rows  of  Eq  (7) 


Yt  " 80  * elXt  * Et 


Yt-1  ■ 60  * SlXt-l  * Et-1 


along  with  p = 1 in  Eq  (8)  yields 


't  • Yt-i  ■ ei  (xt  - xt-i>  * “t 


AYt  = eiAXt  + ut 


.*«*.«•  * »> 


GOR/SM/76D-9 


Eq  (18)  has  independent  error  terms  p^.  Ordinary  Least 
Squares  can  be  performed  on  Eq  (18)  yielding  an  estimate  of 
of  the  original  equation,  Eq  (15).  If  the  assumption  of 
the  first  order  model  is  correct,  and  if  p is  close  to  1, 
then  the  first  difference  estimate  of  B^  will  be  better  than 
the  Ordinary  Least  Squares  estimate.  Also  the  corresponding 
first  differences  prediction  of  Y will  be  better  than  the 
Ordinary  Least  Squares  estimate  of  Y (Ref  8) . 

This  concludes  the  discussion  of  techniques  used  when  a 
first  order  autoregressive  model  is  assumed:  Cochrane-Orcutt , 

Durbin  Nonlinear,  Durbin  Two  Step,  Search,  Bayesian,  and 
First  Differences.  These  techniques  along  with  their  advan- 
tages and  disadvantages  are  summarized  in  Fig.  4. 


Use  Per  Cent  Change  Regression 

Using  per  cent  change  regression  is  not  recommended  in 
textbooks  or  journals  as  one  of  the  procedures  to  use  if 
autocorrelation  is  detected.  However,  it  is  mentioned  here 
because  ASD  used  it  in  their  HOC  report.  See  Chapter  II. 
This  technique  is  similar  to  first  differences,  but  unlike 
first  differences  it  really  is  a complete  respecification  of 
the  model  rather  than  a correction  for  autocorrelation. 
However,  per  cent  changes  may  be  a useful  model  when  the  er- 
rors of  the  original  model  were  not  only  autocorrelated  but 
their  size  also  increased  when  the  independent  and  dependent 
variables  grew  larger. 

This  concludes  the  discussion  of  techniques  that  can  be 
used  when  autocorrelation  is  detected:  ignore  it,  try  a dif- 
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Fig.  4 Techniques  to  use  in  presence  of  aut-  correlation 
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ferent  specification  of  the  model,  assume  a first  order  auto- 
regressive model,  or  use  per  changes.  The  next  chapter  will 
investigate  the  formulation  of  the  historical  indices. 
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IV.  Historical  Indices 
Preface 

One  basic  goal  of  this  thesis  is  to  compare  different 
techniques  that  may  be  used  to  predict  six  aeronautical  price 
indices.  Before  the  techniques  were  compared,  it  was  thought 
that  the  historical  price  indices  should  be  checked  for  va- 
lidity. Invalid  indices  could  effect  the  results  of  the 
thesis.  The  validity  of  the  historical  price  indices  depends 
upon  the  methodology  used  to  calculate  the  indices,  and  the 
accuracy  of  the  calculations.  In  this  thesis  the  methodology 
used  to  calculate  the  historical  price  indices  will  be  as- 
sumed to  be  correct.  However,  the  calculations  will  not  be 
assumed  to  be  correct.  Whenever  feasible,  the  historical  in- 
dices will  be  recalculated.  This  chapter  consists  mainly  of 
examples  of  these  recalculations.  The  recalculated  price  in- 
dices are  listed  in  full  in  Appendix  One.  The  recalculations 
were  done  not  only  to  check  the  validity  of  the  indices,  but 
also  to  gain  a better  understanding  of  the  innerworkings  of 
the  ASD  forecasts. 

Weights  Used 

The  following  weighted  averages  were  used  by  ASD  to  cal- 
culate the  six  composite  aeronautical  price  indices. 


x 
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Airframe  (AF) 


(.4)  AF  Raw  Material  + (.225)  Aircraft  Parts 

§ Equipment 


+ (.375)  7.2%  = Materials 

Overhead  (19) 


(.1425)  Manufacturing  + (.2125)  Engineering 
Labor  Labor 

+ (.1924)  Materials  + (.4526)  7.84  Overhead 

= Development  Airframe 

Index  (20) 


(.2010)  Manufacturing  + (.0650)  Engineering 
Labor  Labor 

+ (.3040)  Materials  + (.4300)  7.8%  Overhead 

= Production  Airframe 

Index  (21) 


Engine 


(.4)  Engine  + (.6)  Aircraft  = Purchased 

Raw  Engine  Labor  Parts 

Material  (22) 


(.2701)  Aircraft  + (.3415)  Engineering  + (.1942)  Engine 

Engine  Labor  Raw 

Labor  Material 

+ (.1942)  Purchased  Parts  = Development 

Engine 

Index  (23) 


(.4184)  Aircraft  + (.0416)  Engineering  + (.2700)  Engine 

Engine  Labor  Raw 

Labor  Material 

+ (.2700)  Purchased  Parts  = Production 

Engine 

Index  (24) 
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Avionics 

(.10)  Avionics  + (.35)  Development  + (.55)  Development 

Raw  Avionics  Overhead 

Material  Labor 

= Development 
Avionics 

Index  (25) 


(.20)  Avionics  + (.30)  Production  + (.5)  Production 

Raw  Avionics  Overhead 

Material  Labor 

= Production 
Avionics 

Index  (26) 

The  weights  and  methodology  used  were  from  ASD  workshops,  a 
Rand  report,  and  a Navy  study  (Ref  4:  Ref  17:  Ref  18).  Note, 
within  each  major  area,  airframe,  engine,  avionics,  the  pro- 
duction and  development  indices  are  closely  related.  Each 
of  the  subindices  listed  in  Eqs  (19)  through  (26)  will  now  be 
discussed  in  length. 

Airframe  Raw  Material  Subindex 

The  airframe  raw  material  subindex  is  a weighted  aver- 
age of  the  Wholesale  Price  Indexes  (base  year  1967)  for  8 
metals.  The  methodology  is  based  on  the  1970  Rand  report 
(Ref  4).  The  Wholesale  Price  Index  (WPI)  for  all  8 metals 
for  years  1958  to  1975  was  not  readily  available.  Three 
sources  were  needed.  Monthly  issues  of  the  Wholesale  Price 
Index  and  Indexes  periodical,  a special  computer  listing  of 
selected  WPI's  from  the  Department  of  Labor,  and  WPI's  listed 
in  the  1970  Rand  report  were  all  used  (Ref  19:  Ref  20:  Ref  4). 
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Data  from  the  three  sources  were  checked  where  it  overlapped 
to  ensure  consistency.  The  following  is  a list  of  the  metals 
used,  the  Bureau  of  Labor  Statistics  codes,  the  weights 
used,  and  an  example  of  the  computations  for  the  1970  Air- 
frame Materials  Subindex: 


Metal 

Code 

Weights 

1970  WPI  (1967=100) 

Finished  Steel 

10- ! 3-02 

.02 

114.4 

Stainless  Steel 

10-13-02-53 

.04 

133.7 

Titanium  Sponge 

10-22-01-56 

.07 

95.5 

Aluminum  Sheets 

10-25-01-01 

.29 

110.7 

Aluminum  Rods 

10-25-01-13 

.11 

93.4 

Aluminum  Extrusions 

10-25-01-17 

.20 

120.6 

Wire  and  Cable 

10-26 

.12 

126.7 

Nuts  and  Bolts 

10-81 

.15 

121.9 

1970  Airframe  Raw  Material  Subindex  114.28 


Aircraft  Parts  and  Equipment  Subindex 

This  index  is  actually  a labor  wage  index.  It  is  the 
average  hourly  earnings  for  aircraft  parts  and  equipment 
production  workers.  Bureau  of  Labor  Statistics  code  3723,9 
converted  into  a 1967  base.  The  index  is  developed  by 
dividing  the  average  hourly  earnings  for  each  year  by  the 
average  hourly  earnings  in  1967  and  then  multiplying  by  100. 
The  data  were  obtained  from  Bureau  of  Labor  Statistics  Em- 
ployment and  Earnings  publication  (Ref  21) . The  following 
is  an  example  calculation  for  the  1970  index: 

Average  hourly  earnings  for  1967,  $3.35 
Average  hourly  earnings  for  1970,  $3.99 


28 

I 

f 


4 


r 


GOR/SM/76D-9 


i 
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1970  Aircraft  Parts  § Equipment  Subindex  = 

x 100  = 119.10  (27) 

Manufacturing  Labor  Subindex 

This  index  is  the  average  hourly  earnings  for  aircraft 
production  workers.  Bureau  of  Labor  Statistics  code  3721,  con- 
verted into  a 1967  year  base.  The  index  is  developed  in  the 
exact  same  way  as  the  Aircraft  Parts  and  Equipment  subindex. 
The  data  is  from  the  exact  same  place  (Ref  21) . The  follow- 
ing is  an  example  calculation: 

Average  hourly  earnings  for  1967,  $3.49 
Average  hourly  earnings  for  1970,  $4.17 
1970  Manufacturing  Labor  Subindex: 

x 100  = 119.48  (28) 

Engineering  Labor  Subindex 

This  index  is  the  average  hourly  earnings  for  engineers 
based  on  year  1967.  The  data  were  supplied  to  ASD  from  six 
large  aerospace  contractors.  The  1974  and  1975  earnings 
were  estimated.  No  example  calculations  are  given  because 
it  might  reveal  proprietary  information. 

7.8%  Overhead  Subindex 

This  index  is  simply  7.8%  compounded  annually  back  and 
forth  from  a 1967  base  of  100.  The  7.8%  was  based  on  over- 
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head  data  provided  to  ASD  from  the  same  large  aerospace  con- 
tractors that  supplied  the  engineering  data.  The  following 
is  an  example  calculation: 

7.8%  Overhead  Subindex  for  1970  = (100) (1  + .078)^ 

= 125.27  (29) 

7.3%  Overhead  Subindex 

This  index  is  simply  7.3%  compounded  back  and  forth 
from  a 1967  base  of  100.  Mr.  Benoy  of  ASD  stated  the  7.3% 
figure  was  based  on  vendors  having  lower  indirect  labor  costs 
than  prime  contractors  (Ref  22) . The  following  is  an  example 
calculation : 


7.3%  Overhead  Subindex  for  1964  = '(100)  (1  + .073) 

= 80.95  (30) 

Engine  Raw  Material  Subindex 


This  index  is  a weighted  average  of  prices  for  8 metals 
converted  into  a 1967  base.  The  prices  were  determined  by 
averaging  weekly  prices  of  the  metals  listed  in  I ron  Age  mag- 
azine. The  metals  selected  and  the  weights  used  were  based 
upon  ASD  analysis  of  data  from  United  States  Air  Force  Air- 
craft Systems  Program  Offices.  Calculation  of  this  index  is 
a very  long  process.  Thus  only  the  1975  figures  were  par- 
tially checked.  The  rest  of  the  data  were  assumed  to  be  cor- 
rect as  listed  in  ASD's  HOC  report  (Ref  3).  The  following 
is  a list  of  the  metals  used,  the  weights  used,  and  the  1975 
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average  prices: 


Metal 

Weights 

1975  Prices 

Titanium 

.329 

225.0 

Chromium 

.165 

235.9 

Aluminum 

.092 

39.6 

Cobalt 

.041 

397.5 

Molydenum 

.030 

259.9 

Vanadium 

.009 

255.1 

Magnesium  Ingot 

.007 

81.8 

Nickel 

.327 

204.5 

1975 

weighted  average 

210.43 

(*/lb) 


Engine  Raw  Material  Subindex  for  1975 

= (1975  prices/1967  prices)  x 100 
= (210.43/105.12)  x 100 
= 200.18  (31) 


Aircraft  Engine  Labor  Subindex 


This  index  is  the  average  hourly  earnings  for  aircraft 
engine  and  engine  parts  production  workers,  Bureau  of  Labor 
Statistics  code  3722.  The  index  is  developed  by  dividing 
the  average  hourly  earnings  for  each  year  by  the  average 
hourly  earnings  in  1967  and  then  multiplying  by  100.  The 
data  were  obtained  from  the  Bureau  of  Labor  Statistic  Em- 
ployment and  Earnings  publication  (Ref  21).  The  following 
is  an  example  calculation: 

Average  hourly  earnings  for  1967,  $3.42 
Average  hourly  earnings  for  1970,  $4.10 
1970  Aircraft  Enginer  Labor  Subindex  = 


x 100  = 119.88 


(32) 
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Avionics  Raw  Material  Subindex 

All  the  Avionics  subindices  are  based  on  methodology 
developed  in  a 1970  United  States  Navy  report  (Ref  17).  This 
index  is  simply  the  arithmetic  average  of  two  Wholesale  Price 
Indexes  (WPI) , base  year  1967.  WPI  1172,  Integrating  and 
Measuring  Instruments,  and  WPI  1178,  Electronic  Components 
and  Accessories.  The  following  is  an  example  calculation: 

1975  Avionics  Raw  Material  Subindex  = 

( . 5)  (1172  WPI)  + (.5)  (1175  WPI)  = .5(140.0)  + 

.5(115.5)  = 127.8  (33) 

Development  Avionics  Labor  Subindex 

This  index  is  a weighted  average  of  hourly  wages  for 
engineers,  engineering  technicians,  draftsmen,  and  produc- 
tion workers,  converted  into  a 1967  base.  ASD  used  the  wage 
data  supplied  or  estimated  in  the  1974  Navy  report,  except 
for  1974  and  1975  where  they  estimated  the  engineers'  wages 
themselves  (Ref  18) . The  following  is  an  example  calcula- 
tion : 


1975  weighted  wages  = .5 (engineer  ave  hrly  wage)  + 

. 2 (engineering  technician)  + . 1 5 (draftsmen)  + 

. 15 (production  workers)  = .5(10.78)  + .2(6.3)  + 
.15(6.1)  + .15(4.66)  = 8.264  (34) 

1975  Development  Avionics  Labor  Subindex  = 

(1975  wages/1976  wages)  * (100)  = 

(8.264/5.0975)  * (100)  = 162.12  (35) 
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Production  Avionics  Labor  Subindex 

This  index  uses  the  same  wages  the  Development  Avionics 
Labor  subinde.- , but  uses  different  weights.  The  following 
is  an  example  calculation: 

1975  weighted  wages  = . 2 (engineers)  = . 1 5 (engineering 
technicians)  + . 05 (draftsmen)  + . 6 (production 
workers)  = .2(10.78)  + .15(6.3)  + .05(6.1)  = 

.6(4.66)  = 6.202  (36) 

1975  Production  Avionics  Labor  Subindex: 

(1975  wages/1967  wages)  * (100)  = 

(6.202/3.804)  x (100)  = 163.04  (37) 

Development  Avionics  Overhead  Subindex 

This  index  is  simply  the  Development  Avionics  Labor  sub- 
index times  1.3274%  compounded  annually  back  and  forth  from 
1967.  Thus  overhead  as  a percentage  of  Avionics  Labor  was 
simply  assumed  to  increase  each  year  by  1.3274%.  The  metho- 
dology is  taken  from  the  1970  Navy  report  (Ref  17).  The  fol- 
lowing is  an  example  calculation: 

1975  Development  Avionics  Overhead  Subindex  = 

162.12  (1  + . 013274) 8 = 180.16  (38) 

Production  Avionics  Overhead  Subindex 

This  index  is  simply  the  Production  Avionics  Labor  sub- 
index times  1.3274%  compounded  annually  back  and  forth  from 
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1 

1967.  Again  overhead  as  a percentage  of  Avionics  Labor  was 
assumed  to  increase  each  year  by  1.3274%.  The  methodology 
is  again  from  the  Navy  study  (Ref  17) . The  following  is  an 
example  calculation: 

1975  Production  Avionics  Overhead  Subindex  = 

163.04  (1  + . 013274) 8 = 181.18  (39) 

That  concludes  a brief  explanation  of  the  methods  used 
in  the  preparation  of  the  historical  values.  The  recalcu- 
lated indices  are  listed  in  Appendix  A.  The  recalculated  in- 
dices compare  very  closely  to  the  ASD  calculations  in  the 
HOC  report  with  one  exception  (Ref  4)  . The  ASD  Airframe 
Raw  Material  Index  value  was  slightly  lower  than  the  recal- 
culated values  prior  to  1967.  This  was  due  to  differences 
in  the  ASD  calculations  when  the  Bureau  of  Labor  Statistics 
switched  from  a 1957-1958  base  to  a 1967  base.  ASD  has  re- 
calculated the  Airframe  Raw  Material  Subindex  values  for 
use  in  future  reports. 

One  can  see  that  there  are  basically  three  sets  of  in- 
dices: labor  wage  indices,  raw  material  price  indices,  and 

overhead  cost  indices.  The  labor  wage  indices  were  based  on 
average  hourly  earnings.  The  raw  material  indices  were  a 
combination  of  Wholesale  Price  Indexes,  and  the  overhead  in- 
dices increased  at  a given  percentage  each  year.  Many  weight- 
ed averages  were  used  to  construct  the  historical  indices. 

It  is  clear  that  one  could  spend  the  entire  time  in  a thesis 
discussing  the  choice  and  development  of  price  indices.  In  - * 
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fact,  a 1971  Air  Force  Institute  of  Technology  thesis  did 
just  that  (Ref  23).  However,  it  was  decided  not  to  inves- 
tigate the  historical  indices  any  further.  The  rest  of  the 
thesis  will  be  spent  investigating  the  prediction  of  the  con- 
structed indices.  First,  however,  the  independent/predictor 
variables  will  be  discussed. 

Independent/Predictor  Variables 

There  are  five  variables  in  this  thesis  that  were  used 
or  tried  as  independent/predictor  variables  in  the  regres- 
sions. These  variables  are:  GNP  price  deflator,  Primary 
Metals  Index,  Electrical  Machinery  Index,  the  money  supply 
Ml,  and  the  U.S.  total  unemployment  rate.  The  aircraft  in- 
dustry unemployment  rate  would  have  been  used  instead  of  the 
total  U.S.  unemployment  rate,  but  data  were  only  available 
on  the  aircraft  unemployment  rate  after  1964  (Ref  24).  His- 
torical values  for  the  Primary  Metals  Index  and  the  Electric 
Machinery  Index  were  taken  from  the  ASD  HOC  report  (Ref  4) . 

The  original  source  was  Wharton  Econometric  Forecasting  As- 
sociates, Inc.  Historical  values  for  the  GNP  price  deflator. 

Ml,  and  the  national  unemployment  rate  were  taken  from  the 
January  1976  Economic  Report  of  the  President  (Ref  25) . The 
historical  independent  variables  are  all  listed  in  Appendix 
B.  Predictions  for  some  of  these  variables  are  also  listed 
in  Appendix  B.  The  predictions  are  from  the  Wharton  Econo- 
metric Forecasting  Associates,  Inc.  More  will  be  said  about 
these  predictions  in  Chapter  VII. 
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This  concludes  a discussion  of  the  historical  data. 

The  next  chapter  will  compare  a technique  used  in  the  HOC 
report,  subindex  per  cent  change  regression,  with  a similar 
technique,  aggregate  per  cent  change  regression. 
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V.  Per  Cent  Change  Regressions 
Preface 

This  chapter  describes  two  distinct  approaches  to  index 
prediction  using  per  cent  change  regression.  The  first  ap- 
proach is  the  subindex  method  used  by  ASD  in  the  HOC  report. 

See  Fig.  1.  The  regression  equations  using  this  approach 

2 

are  shown  below  along  with  R , the  coefficient  of  determina- 
tion, the  t statistics  for  the  coefficients,  DW,  the  Durbin 

2 

Watson  statistic  for  autocorrelation,  and  S , the  calculated 
variance.  Predictions  and  prediction  intervals  will  also 
be  shown. 

The  second  approach  involves  first  aggregating  the  data, 

and  then  performing  the  regressions.  The  regression  equa- 

2 

tions  for  this  second  approach  will  be  shown  along  with  R , 

2 

the  t statistics,  DW,  and  S . The  two  approaches  are  com- 
pared in  the  end  of  the  chapter.  All  of  the  calculations  in 
this  chapter  were  performed  on  the  indices  recalculated  in 
Chapter  IV. 

Subindex  Method 


Regression  Equations 

A % Engineering  Labor  = 1.582  + .8286  A % GNP, 

(1.76)  (4.27)  a 

DW  = 2.01  R2  = .55  S2  = 3.8739  (40) 
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A % Manufacturing  Labor  = 2.124+  .8602  A % GNP, 

(2.99)  (5.60)  a 


DW  = 2.40 


R = .68 


= 2.4342 


A I AC  Engine  Labor  = 2.281  + .7748  A % GNP, 

(3.63)  (5.71)  a 


DW  = 2.15 


R = .68 


= 1.8784 


A % AC  Parts  8 Equipment  = 2.321  + .6673  A % GNP, 

(5.40)  (7.19)  a 


DW  = 2.37 


R = .77 


S‘  = .8888 


A % Development  Avionics  Labor  = 3.445  + .4429  A % GNP, 

(6.40)  (3.81)  Q 


DW  = 1.78 


R = .49 


= 1.3955 


A % Avionics  Production  Labor  = 2.914  + .5361  A I GNP, 

(5.35)  (4.55)  Q 

DW  = 2.44  R2  = .58  S2  = 1.4281 


A 


AF  Raw  Material  = -1.381  + 1.143  A \ P Metal 

(-.87)  (4.46) 

DW  = 1.98  R2  = .57  S2  = 24.461 


A 


Engine  Raw  Material 
DW  = 1.82  R2  = 


-1.516  + 1.296  A \ P Metal 
(-.96)  (5.05) 

63  S2  = 24.565 


A % Avionics  Raw  Material  = .7964  + .6907  A % E Mach 

(1.72)  (4.93) 

DW  = 2.25  R2  = .62  S2  « 3.5308 


where  A % 


is 


Indext  - Indext  ^ 
Indext  ^ 


x 100 


(41) 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


(48) 
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GNP,  is  GNP  Price  Deflator 
a 

PMetal  is  Primary  Metals  Index 

E Mach  is  Electrical  Machinery  Index 

Numbers  in  parenthesis  are  t statistics  for  the  coef- 
ficients 

All  equations  except  that  for  Airframe  Raw  Material  are  very 
similar  to  the  equations  listed  in  ASD  report  HOC  (Ref  3: 
C-l).  The  difference  in  the  Airframe  Raw  Material  equations 
was  a direct  result  of  differences  in  calculation  of  the  his- 
torical data.  See  Chapter  IV.  The  corresponding  equation 
listed  in  report  110C  is: 

A % AF  Raw  Material  = 1.289  + 1.146  A % PMetal  (49) 

For  all  regression  equations,  visual  residual  analysis 
did  not  indicate  any  strong  departure  from  the  assumptions 
of  Ordinary  Least  Squares.  The  Durbin  Watson  statistics  all 
tend  to  indicate  that  the  errors  were  not  autocorrelated . 
However,  four  coefficients  were  found  not  to  be  significant- 
ly different  from  0 at  the  95%  confidence  level.  The  four 
coefficients  were  the  constants  in  each  of  the  raw  material 
regressions  and  the  constant  in  the  engineering  labor  regres- 
sion. The  regression  equations  will  be  used,  as  is,  with 
the  constant  terms  present,  to  predict  future  values.  This 
is  done  so  that  a direct  comparison  can  be  made  with  the  110C 
methodology 

Predictions 


Changing  per  cent  change  predictions  into  subindex  pre- 
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dictions  is  a simple  process.  The  per  cent  change  predic- 
tions are  repeatedly  applied  to  the  last  known  index,  the 
1975  value.  The  following  is  an  example  of  the  equations 
used : 


A % 7 s 

76  index  = 7 5 index  (1  + ) 

A % 7,  A % 77 

77  index  = 75  index  (1  + 1QJ^  ) (1  + y^7  - -) 


(50) 

(51) 


where  A%y^  stands  for  the  per  cent  change  prediction  for 
1976.  The  predictions  for  the  subindices  are  listed  in  Table 
I.  The  predictions  for  all  subindices  are  almost  identical 
to  their  counterparts  in  the  HOC  report  (Ref  3:12-14).  Af- 
ter the  subindices  are  predicted,  they  are  weighted  and  ag- 
gregated using  Eqs  (19)  through  (25) . This  yields  the  pre- 
dictions for  the  final  six  indices.  These  predictions  are 
listed  in  Table  II.  No  direct  comparison  of  these  final  six 
indices  can  be  made  with  the  110C  report.  The  final  six  in- 
dices are  listed  in  the  HOC  report  only  in  a Fiscal  Year 
base . 


Prediction  Intervals 

Only  95%  prediction  intervals  will  be  given  in  this 
thesis.  If  the  regression  model  is  correct  and  continues  to 
hold  into  the  future,  and  if  the  Wharton  predictors  are  known 
without  error,  then  the  following  could  be  said  about  the 
prediction  intervals.  If  one  could  somehow  collect  new  his- 
torical data  on  the  same  indices,  do  the  regressions  again, 
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Table  I 


Subindex  Predictions  and  Prediction  Intervals 
Using  Per  Cent  Change  Regressions 


CY 

Engineering 

Manufacturing 

AC  Parts  5 

Labor 

Labor 

Equipment 

1976 

182. 97±  7.09 

190 . 48±  5 . 81 

175.12±3.25 

1977 

195.17+10.39 

204 . 59±  7 . 00 

186.36±4.86 

1978 

210.10 

221.81 

199.79 

1979 

227.09 

241.50 

214.89 

1980 

242.98 

260.21 

229.25 

1981 

256.51 

276.50 

241.93 

1982 

273.39 

296.72 

257.29 

1983 

290.94 

317.92 

273.28 

1984 

304.68 

335.03 

286.53 

1985 

316.20 

349.81 

298.26 

CY 

AC  Engine 
Labor 

Development 
Avionics  Labor 

Production 
Avionics  Labor 

1976 

188.44+5.09 

171. 9 6± 4. 01 

172.97i4.-08 

1977 

201.70+7.48 

182 . 56±  5.87 

183.71±5.97 

1978 

217.74 

194.77 

196.27 

1979 

235.95 

208. 26 

210.25 

1980 

253.28 

221.46 

223.74 

1981 

268.49 

233.81 

236.03 

1982 

287.16 

248.11 

250.54 

1983 

306.69 

263.07 

265.68 

1984 

322.69 

276.55 

278.82 

1985 

336.69 

289.33 

290.92 

CY 

AF  Raw 

Engine  Raw 

Avionics  Raw 

Material 

Material 

Material 

1976 

175.62117.00 

213.55120.55 

135.9015.62 

1977 

188. 31±24. 91 

231.15130.26 

143.8918.12 

1978 

201.77 

250.00 

151.78 

1979 

215.04 

268.76 

159.60 

1980 

228.80 

288.40 

167.35 

1981 

239.66 

304.07 

175.06 

1982 

254.56 

325.66 

182.73 

1983 

269.38 

347.32 

190.38 

1984 

278.67 

361.06 

198.00 

1985 

284.66 

370.04 

205.62 

41 


GOR/SM/76D-9 


Table  II 


Composite  Index  Predictions  and  Prediction  Intervals: 
Using  Subindex  Per  Cent  Change  Regressions 


CY 

Development 

Airframe 

Production 

Airframe 

1976 

189 . 70±1 . 92 

189.54+1.83 

1977 

203 . 7 0±  2 . 81 

203.64±2.69 

1978 

219.49 

219.43 

1979 

236.79 

236.60 

1980 

254.44 

254.25 

1981 

271.71 

271.58 

1982 

291.50 

291.43 

1983 

312.49 

312.47 

1984 

332.50 

332.56 

1985 

352.43 

352.59 

CY 

Development 

Engine 

Production 

Engine 

1976 

193.40±5.89 

197. 70±  7.13 

1977 

207. 4 8± 8. 66 

212.56+10.50 

1978 

223.90 

229.62 

1979 

241.85 

247.98 

1980 

259,31 

266.13 

1981 

274.07 

281.44 

1982 

292.93 

301.14 

1983 

312.36 

321.39 

1984 

326.97 

336.45 

1985 

338.76 

348.44 

CY 

Development 

Avionics 

Production 

Avionics 

1976 

180.27+3.79 

176.46±4.00 

1977 

192 . 84±5 . 59 

188 . 69±5 . 87 

1978 

207.19 

202.69 

1979 

223.03 

218.14 

1980 

238.83 

233.38 

1981 

247.39 

247.77 

1982 

271.42 

264.38 

1983 

289.79 

281.83 

1984 

306.92 

297.69 

1985 

323.59 

312.83 

and  each  time  calculate  95%  predictions  intervals,  then  the 
probability  is  that  95  out  of  100  of  the  prediction  inter- 
vals calculated  would  contain  the  actual  future  value  of  the 
index  (Ref  25:304).  However,  one  cannot  collect  the  histori- 
cal data  again.  In  practice,  then,  the  following  is  said: 
There  is  a 95%  probability  that  the  actual  future  index 
values  will  fall  within  the  prediction  intervals.  This  of 
course,  assumes  that  the  relationship  between  the  predictor 
variables  and  the  indices  will  continue  to  be  the  same,  and 
that  the  Wharton  Predictions  are  without  error.  The  effect 
of  errors  in  the  Wharton  predictors  will  be  discussed  in 
Chapter  VII. 

Calculation  of  prediction  intervals  is  a difficult  pro- 
cess. In  fact  the  calculations  required  estimates  of  the 
covariances  between  the  subindices  and  the  assumption  of  a 
near  normal  distribution  for  the  predictions.  Since  the  cal- 
culations were  burdensome,  only  prediction  intervals  for  1976 
and  1977  were  calculated.  The  prediction  intervals  are  listed 
in  Tables  I and  II.  Prediction  intervals  for  future  years 
will  grow  progressively  larger  and  larger  at  a rate  compar- 
able to  the  aggregate  method  (See  Table  III).  The  full  ex- 
planation of  predicion  interval  construction  is  contained  in 
Appendix  C. 
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Aggregate  Method 


Regression  Equations 

A % Air  Frame  Development  = 3.959  + .5544  A \ GNP, 

(13.98)  (9.06)  a 


DW  = 2.36 


R = .845 


= .3858 


A % Air  Frame  Production  = 3.753  + .6020  A % GNP 

(15.60)  (11.54)  c 


DW  = 2.80 


R = .900 


= .2785 


A % Engine  Development  = -.9228  + 1.4054  A % GNP 

(-1.70)  (11.97)  c 


DW  = 2.42 


1.4210 


A % Engine  Production  = -1.825  + 1.6207  A % GNP, 

(-3.44)  (14.13)  a 


DW  = 2.17 


R = .930 


S = 1.3549 


A % Avionics  Development  = 3.3155  + .5521  A % GNP, 

(6.52)  (5.03)  a 


DW  = 1.74 


R = .627 


S'  = 1.2435 


A % Avionics  Production  = 2.0740  + .7087  A \ GNP, 

(4.37)  (6.71)  a 


DW  = 2.28 


R“  = .761 


= 1.0850 


where  all  symbols  have  the  same  meaning  as  used  in  the  sub- 
index method 


For  all  regression  equations,  visual  residual  analysis  did 
not  indicate  any  strong  departure  from  the  assumptions  of 
Ordinary  Least  Squares.  All  the  Durbin  Watson  statistics, 
except  for  Production  Airframe,  tended  to  indicate  that  the 
errors  were  not  autocorrelated . The  Durbin  Watson  statistic 
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for  the  Airframe  regression  was  in  the  inconclusive  region. 
Only  one  coefficient,  the  constant  term  for  the  Development 
Engine  equation,  was  not  significantly  different  from  zero  at 
the  95°s  confidence  level.  ASD  did  not  use  the  aggregate  per 
cent  change  method  in  any  of  their  reports. 

Predictions 

The  predictions  were  done  using  Eqs  (50)  and  (51)  . No 
other  steps  are  needed  since  the  data  has  already  been  aggre- 
gated. The  predictions  are  listed  in  Table  III. 

Prediction  Intervals 

Unlike  the  subindex  method,  the  calculation  of  the  pre- 
diction intervals  did  not  require  estimates  of  covariances; 
however,  a near  normal  distribution  was  assumed  for  the  pre- 
dictions. The  prediction  intervals  for  all  six  composite 
indices  are  listed  in  Table  III.  The  prediction  interval  for 
Production  Avionics  is  shown  in  Fig.  5.  The  explanation  of 
the  prediction  interval  construction  for  the  aggregate  method 
is  shown  in  Appendix  C. 

Comparison  of  the  Aggregate  and  the  Subindex  Methods 

It  is  important  to  realize  that  the  difference  between 
the  subindex  per  cent  change  regression  and  the  aggregate  per 
cent  regression  method  is  not  just  the  aggregation.  There 
are  differences  in  the  independent  variables  and  in  the  treat- 
ment of  overhead  between  the  two  methods.  The  subindex  method 
uses  the  Primary  Metals  Index  as  the  independent  variable  for 


45 


GOR/ SM/76D- 9 


Table  III 


Composite  Index 
Using  Aggregate 

Predictions  and 
Per  Cent  Change 

Prediction  Intervals: 
Regressions 

CY 

Development 

Airframe 

Production 

Airframe 

1976 

189. 77± 

2.30 

189 .71  + 1.95 

1977 

203. 74± 

3.39 

203 . 85±  2 . 88 

1978 

220. 08± 

4.37 

220.48±3.71 

1979 

238. 37± 

5.31 

239 . 18±4 . 52 

1980 

256. 45± 

6.18 

257 . 57±  5 . 26 

1981 

273. 43± 

7.02 

274.69i5.98 

1982 

293. 41± 

7.89 

294 . 98±  6 . 7 2 

1983 

314. 52± 

8.78 

316.41+7.48 

1984 

333. 58± 

9.67 

335.51+8.25 

1985 

351. 70± 

10.60 

353 . 46±9 . 04 

CY 

Development 

Production 

Engine 

Engine 

. 1976 

194. 70±  4.53 

199. 56±  4.52 

1977 

209. 71±  6.67 

215. 78±  6.66 

1978 

22-9. 34±  8.61 

235. 43±  8.63 

1979 

252.53+10.54 

263 . 31±10 . 60 

1980 

273.40+12.35 

286.41±12.47 

1981 

289.35+14.09 

303.48±14.26 

1982 

311 . 21± 1 5 . 88 

327.62+16.12 

1983 

333 . 86±17 . 71 

352. 6 2± 18. 03 

1984 

348 . 56±19 . 53 

307.85+19.93 

1985 

358 . 36±  21.35 

376. 97±21 . 82 

CY 

Development 

Production 

Avionics 

Avionics 

1976 

180. 49±  3.96 

175.42+  3.60 

1977 

192. 59±  5.80 

186.69+  5.27 

1978 

206. 76±  7.45 

200.25+  6.76 

1979 

222.58+  9.11 

215.54+  8.18 

1980 

237 . 99±1 0 . 54 

230. 00±  9.49 

1981 

252. 20±11. 90 

242 . 60±10 . 71 

1982 

268.96+13.31 

258 . 01±11 . 98 

1983 

283 . 85±14 . 73 

274.03+13.25 

1984 

299 . 20±16 . 13 

287. 08±14 . 51 

1985 

313. 51±17. 57 

298.43+15.79 
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the  airframe  and  engine  material  regressions,  and  the  Elec- 
tric Machinery  Index  as  the  independent  variable  in  the  avi- 
onics materials  regression.  The  GNP  price  deflator  is  used 
in  the  other  subindex  regressions.  In  the  aggregate  method, 
only  the  GNP  price  deflator  is  used  as  the  independent  vari- 
able in  the  regressions.  Primary  Metals  and  Electrical  Mach- 
inery were  tried  as  additional  independent  variables  in  the 
aggregate  method,  but  they  were  found  not  to  be  significant 
at  the  951  confidence  level. 

In  the  subindex  method,  overhead  is  not  regressed  but 
rather  projected  at  a predetermined  rate.  For  the  airframe 
indices  two  different  overhead  subindices  were  used.  One 
overhead  subindex  is  7.3%  of  the  value  for  the  previous  year. 

The  other  overhead  subindex  is  7.8%  of  the  value  of  the  pre- 
vious year.  The  avionics  overhead  as  a percentage  of  avionics 
labor  increases  each  year  by  1.3274%.  In  the  subindex  method 
these  overheads  are  projected  into  the  future  at  their  re- 
spective rates.  In  the  aggregate  method  the  overhead  sub- 
indices are  first  aggregated  with  the  other  subindices  and 
then  regressed,  and  then  predicted.  See  Fig.  1. 

Keeping  the  procedural  differences  in  mind,  the  predic- 
tions for  both  the  subindex  and  aggregate  method  can  be  ana- 
lyzed. The  predictions  of  the  airframe  indices  differed 
little  between  the  methods.  The  aggregate  method  engine  in- 
dices had  slightly  higher  values  than  their  subindex  counter- 
parts. The  predictions  for  the  avionics  indices  were  signif- 
icantly lower  using  the  aggregate  method  than  the  subindex 
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method.  The  reason  for  the  differences  in  the  engine  indices 
is  not  known.  Part  of  the  difference  in  the  avionics  in- 
dices is  probably  due  to  the  dramatic  difference  in  overhead 
calculation  between  the  two  methods. 

Both  methods  have  their  own  advantages  and  disadvantages 
The  aggregate  method  is  computationally  simpler  and  quicker. 
Only  six  regressions  are  done  instead  of  nine.  The  weights 
used  in  report  HOC  are  only  estimates  at  best.  The  subin- 
dex method  would  permit  a user  of  the  report  to  use  his  own 
weights  to  develop  production  and  development  composite  in- 
dices without  having  to  do  any  further  regressions.  In  the 
aggregate  method  a user  who  wants  to  use  his  own  weights  has 
to  do  the  regressions  over.  However,  as  noted  below,  this 
probably  presents  no  significant  disadvantage.  The  predic- 
tion intervals  give  little  advantage  to  one  method  over  the 
other.  Three  indices  have  smaller  prediction  intervals  for 
the  aggregate  method,  and  three  indices  have  smaller  predic- 
tion intervals  for  the  subindex  method. 

Given  three  possible  independent  variables,  the  GNP 
price  deflator,  the  Primary  Metals  Index,  and  the  Electrical 
Machinery  Index,  and  given  the  choice  between  the  per  cent 
change  subindex  method  and  the  aggregate  per  cent  change  meth 
od,  the  author  recommends  the  aggregate  method.  Computation- 
ally it  is  simpler  and  quicker.  In  addition,  if  a user  were 
inclined  to  use  his  own  weights  he  would  probably  want  to 
make  other  changes,  and  might  even  do  an  entire  report  him- 
self. Furthermore,  the  author  prefers  the  treatment  of  over- 
head used  in  the  aggregate  method. 
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VI.  Improving  the  Estimates 

This  chapter  investigates  different  techniques  that 
might  be  used  to  imporove  ASD's  estimation  of  the  six  com- 
posite aeronautical  price  indices.  Techniques  including 
simple  regression,  regression  with  nonlinear  variables,  mul- 
tivariate regression.  Generalized  Least  Squares,  First  Dif- 
ferences, and  Per  Cent  Changes  will  be  discussed. 

Defining  the  Beginning 

Before  the  investigation  began,  it  was  necessary  to  de- 
fine a starting  point.  The  historical  subindices  could  have 
been  used  as  the  starting  point  or  the  subindex  data  could 
have  been  aggregated  and  the  composite  indices  could  have 
been  used  as  the  starting  point.-  Unfortunately,  th.ere  was 
not  enough  time  to  do  parallel  investigations,  one  for  each 
possible  starting  point.  The  composite  indices  were  chosen 
as  the  starting  point  for  three  reasons:  first,  fewer  re- 

gressions were  needed;  second,  ASD  was  contemplating  using 
the  composite  technique  in  their  next  report;  and  third,  the 
composite  technique  was  recommended  in  Chapter  V when  dealing 
with  per  cent  change  regressions.  It  is  important  to  realize 
that  the  results  of  this  chapter  are  based  upon  starting  with 
the  composite  data.  Different  relations  that  might  result 
are  discussed  in  the  end  of  this  chapter. 

Simple  Regression 

Regressing  the  six  composite  indices  against  the  GNP 
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price  deflator  is  the  easiest  technique  to  use.  The  results 
are  listed  below. 


Development  Air  Frame  = -54.789  + 1.900  GNP, 

(-10.22)  (30.60)  d 


DW  = .25 


R = .98 


S = 21.02 


Production  AF  = -53.743  + 1.884  GNP 


(-10.67)  (32.30) 


DW  = .24 


R = .98 


S = 18.54 


Development  Engine  = -20.532  + 1.558  GNP 

(-8.24)  (153.93)  C 


DW  = .65 


R = 


S = 4.55 


Production  Engine  = -18.363  + 1.542  GNP 

(-4.27)  (30.91)  C 


DW  = .42 


S4  = 13.56 


Development  Avionics  = -37.050  + 1.684  GNP, 

(-7.80)  (28.46)  a 


DW  = .28 


R = .98 


S4  = 18.85 


Production  Avionics  = -22.855  + 1.521  GNP 

(-7.05)  (40.49)  c 


DW  = .39 


R = .99 


S4  = 7.69 


This  was  basically  the  same  technique  as  used  by  ASD  in 
the  110A  and  HOB  reports.  See  Chapter  II.  The  result  is 
that  the  Durbin  Watson  statistics  all  indicate  that  the  data 
may  be  highly  autocorrelated . This  is  a warning  that  the 
estimates  obtained  using  this  technique  may  not  be  the  best 


available . 
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To  get  a better  idea  of  what  was  actually  happening, 
the  historical  data  were  plotted  along  with  the  regression 
lines  in  Figs.  6 to  11.  Careful  examination  of  these  plots 
yields  interesting  results.  There  is  a definite  pattern  to 
the  data.  In  Figs.  6,  7,  10,  11  the  data  starts  below  the 
regression  line,  remains  there  for  awhile,  and  then  goes 
above  the  line,  and  finally  ends  up  below  the  line  again. 

In  Figs.  8 and  9 the  exact  opposite  happens.  From  these  pat- 
terns one  can  immediately  tell  that  the  Durbin  Watson  statis- 
tic will  be  very  low  indicating  possible  strong  autocorre- 
lation. The  data  would  have  to  be  displaced  randomly  above 
and  below  the  line  for  the  Durbin  Watson  statistic  to  indi- 
cate that  the  data  were  not  autocorrelated . Thus,  even 
though  R is  very  high  and  the  data  points  are  close  to  the 
regression  line,  the  low  Durbin  Watson  statistic  indicates 
that  procedures  other  than  simple  regression  may  be  more  ac- 
curate. Since  accuracy  is  so  important  in  the  price  index 
forecasting,  other  regression  techniques  will  be  investi- 
I gated. 

Regression  with  Nonlinear  Variables 

Would  a nonlinear  model  fit  the  data  better?  A quad- 
ratic model  would  certainly  fit  well  over  a few  years,  but 
then  the  data  swings  back  the  other  way  and  the  quadratic 

[ model  would  no  longer  fit  well.  High  order  curves  may  fit 

\ 

the  data  well,  but  it  would  be  hard  to  justify  using  a high- 
er order  curve.  First,  and  most  important  of  all,  a linear 
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Fig.  6.  Development  Airframe  Index  vs.  GNP  Price  Deflator 
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Production  Airframe  Index  (1967  = 100) 
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Development  Engine  Index  (1967  = 100) 


Production  Engine  Index  (1967  = 100) 
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relationship  between  the  GNP  price  deflator  and  the  composite 
indices  makes  economic  sense.  Higher  order  curves  do  not 
make  economic  sense.  Second,  there  were  not  enough  data  to 
reasonably  estimate  the  parameters  for  the  higher  order 
curve.  So,  no  attempts  were  made  at  fitting  higher  order 
curves  to  the  data. 

Multivariate  Regression 

If  a nonlinear  model  did  not  seem  reasonable,  what 
about  additional  variables?  Three  additional  variables  were 
tried:  the  money  supply  Ml,  lagged  one  year,  the  metals  in- 
dices discussed  in  Chapter  IV,  and  the  national  unemployment 
rate.  The  unemployment  rate  for  the  aircraft  industry  would 
have  been  used,  but  only  data  after  1964  was  available.  The 
results  of  the  regressions,  are  listed  in  Table  IV.  The 
Durbin  Watson  statistics,  the  sign  of  coefficient  for  the 
additional  variable,  and  whether  the  coefficient  was  signifi- 
cant at  the  95%  confidence  level  are  all  listed.  As  expected 
with  the  addition  of  any  variable,  the  Durbin  Watson  sta- 
tistic increases.  In  fact,  in  two  cases  the  Durbin  Watson 
statistic  moved  out  of  the  critical  region  for  acceptance 
of  autocorrelation.  One  case  was  for  Production  Avionics 
Index  with  GNPd  and  the  Electrical  Machinery  Index.  The  sign 
of  the  coefficient  in  this  case  makes  little  economic  sense. 
The  negative  sign  says  when  the  Electrical  Index  goes  up  the 
Production  Avionics  Index  goes  down.  The  other  case  where 
the  Durbin  Watson  statistic  moved  out  of  the  critical  region 
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Table  IV 

Durbin  Watson  Statistics 


GNP  , 

GNP , 

GNP, 

Plu§ 

Plus 

PluS 

Straight  Regressions  GNP^  Metal 

Ml 

t-1 

Unemployment 

Development 

AF 

.25 

.83$ 

.78? 

1.23$ 

Production 

AF 

. 24 

.80$ 

.71? 

1.25$ 

Development 

Engine 

.65 

1.19$ 

1 

.36$ 

1.63+  j 

Production 

Engine 

.42 

.90+ 

1 

.12$ 

s 

1.06  + 

Development 

Avionics 

.28 

1.41$ 

.85? 

1.15$ 

Production 

Avionics 

.39 

1.62- 

1 

.08  + 

1.36$ 

Du 

1.39 

1.53 

1 

.53 

1.53 

Du  is  from  the  Durbin  Watson  tables 

+ , - indicate  the  sign  of  the  coefficient  for  the  ad- 
ditional variable 

s means  the  coefficient  of  the  additional  variable  was 
statistically  significant  at  the  95%  confidence  level. 

N means  the  coefficient  of  the  additional  variable  was 
not  statistically  significant. 


was  for  Development  Engine  Index  with  GNP^  and  unemployment. 
One  might  be  able  to  justify,  economically,  the  addition  of 
unemployment  if  one  ignored  the  signs  of  the  coefficients  of 
the  other  indices.  The  signs  indicate  that  unemployment  is 
positively  related  to  the  Engine  indices  and  negatively  re- 
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lated  to  the  Airframe  and  Avionics  indices.  It  would  be 
hard  to  explain  this  phenomenon  in  economic  or  business 
terms . 

The  difference  in  signs  of  the  coefficients  is  not  hard 
to  explain,  statistically.  An  additional  variable  will  be 
statistically  significant  if  it  relates  well  to  the  resid- 
uals, the  differences  between  the  regression  line  and  the 
actual  values,  of  the  GNP^  only  regress  ion . As  shown  in 
Figs.  6 to  11,  the  pattern  of  the  residuals  for  the  Airframe 
and  Avionics  indices  is  opposite  that  of  the  Engine  indices. 
The  sign  of  the  coefficients  reflect  these  same  patterns. 

One  can  even  predict  whether  or  not  an  additional  vari- 
able will  be  statistically  significant.  The  pattern  of  the 
residuals  shows  that  any  additional  variable  that  has  a dip 
or  rise  in  the  middle  of  its  data  will  probably  be  statis- 
tically significant.  This  is  why  one  must  be  so  careful 
when  doing  a regression.  Many  variables  could  be  added  to 
a regression  that  are  statistically  significant,  but  which 
are  not  consistent  with  the  underlying  phenomena  being  mod- 
eled. None  of  the  variables  tried  here  are  recommended  as 
additional  variables  in  Ordinary  Least  Squares  analysis  on 
the  composite  indices. 

First  Order  Autoregressive  Model 

The  next  technique  tried  was  the  assumption  of  a first 
order  autoregressive  model.  The  Durbin  Two  Step,  the  Search, 
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and  the  Bayesian  techniques  were  used  to  try  and  estimate  p. 
These  techniques  are  explained  in  Chapter  III.  Computer  code 
for  the  Search  and  Bayesian  techniques  is  included  in  Appen- 
dix D.  The  results  of  the  different  techniques  are  dis- 
played in  Table  V.  As  explained  in  Chapter  III,  the  Search 
answers  were  considered  the  most  reliable,  the  Durbin  Two 
Step  answers  were  used  as  a check,  and  the  Bayesian  answers 
were  used  to  estimate  how  much  greater  than  one  p might  be. 
There  was  a problem  with  the  Bayesian  technique.  For  some 
cases  the  variance  did  not  reach  a single  minimum,  rather 
the  variance  reached  a local  minimum  then  increased  and  then 
decreased  to  another  local  minimum.  In  these  cases  only 
the  first  local  minimum  reached  was  displayed.  These  cases 
are  indicated  by  a ^ in  Table  V. 

The  results  of  the  one  variable.  Search,  Durbin  Two 
Step,  and  Bayesian  techniques  are  consistent  with  the  raw 
data.  The  simple  regressions  done  in  the  beginning  of  the 
chapter,  all  resulted  in  very  small  Durbin  Watson  statistics. 
These  small  Durbin  Watson  statistics  are  consistent  with  a 
large  value  of  p in  the  first  order  autoregressive  model. 

The  Development  Engine  index  had  the  largest  Durbin  Watson 
statistic  for  the  simple  regression.  As  expected,  this  in- 
dex had  the  smallest  estimate  of  p for  the  one  variable 
Search  technique.  Also,  a close  look  at  Figs.  6 through  11 
reveals  a slight  tendency  for  the  data  to  be  close  to  the  re- 
gression line  for  the  small  values  of  the  GNP  price  deflator 
and  any  farther  away  for  large  values  of  the  GNP  price  defla- 
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Table  V 
Estimates  of  p 


Two  Step 
Durbin 
GNP^  only 

Search 
GNP^  only 

Bayesian  Search 
GNP,  6 
GNP^  only  Metals 

Search 

GNP  .^General/ 
Unemployment 

Development 

Airframe 

1.04 

>.99 

1.1* 

>.99- 

>.99- 

Production 

Airframe 

1.G3 

>.99 

1.1* 

>.99 

>.99- 

Development 

Engine 

1 . 01 

.95 

.71  + 

.3  + 

Production 

Engine 

1.31 

>.99 

1.1* 

>.99+ 

> .99+ 

Development 

Avionics 

.94 

>.99 

1.1 

.35- 

>.9  - 

Production 

Avionics 

.85 

>.99 

1.1 

.21- 

>.9  - 

* means  there  was  more  than  one  local  minimum  of 
the  variance.  The  first  minimum  reached  is 
listed. 

+ , - are  the  signs  of  the  coefficients  of  the  ad- 
ditional variables. 


tor.  This  apparent  divergence  of  the  data  is  consistent 
with  a p slightly  greater  than  one  in  the  first  order  auto- 
regressive model. 

The  Search  technique  was  also  used  with  an  additional 
variable.  Both  the  general  unemployment  rate  and  the  metals 
indices  were  tried.  Ml  was  not  tried  because  of  its  high 
correlation  with  GNP^  creating  dependency  problems.  See  OLS 
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assumption  4,  Cahpter  III.  The  results  are  displayed  in 
Table  V.  The  signs  of  the  coefficients  of  the  additional 
variables  are  also  listed  in  Table  V.  The  signs  follow  the 
same  pattern  as  they  did  earlier  in  the  ordinary  linear  model 
regressions . 

Partially  because  of  the  signs  of  the  coefficients, 
predictions  were  only  made  for  two  cases.  Predictions  were 
made  for  the  Development  Engine  index  for  p = .95,  and  for  p 
= .71  with  the  Primary  Metals  index  as  an  additional  variable. 
The  predictions  along  with  the  prediction  intervals  are  listed 
in  Table  VI.  Predictions  were  not  made  for  the  Avionics  in- 
dices with  the  addition  of  the  Electrical  Machinery  index. 

This  was  because  the  coefficient  of  the  Electrical  Metals  in- 
dex was  negative  and  made  no  economic  sense.  Predictions 
were  not  made  for  the  Development  Engine  index  with  unemploy- 
ment added.  As  explained  earlier  in  the  chapter,  the  pat- 
tern  of  the  signs  of  the  coefficients  of  unemployment  in  the 
six  indices  is  difficult  to  explain  economically.  The  pre- 
dictions were  only  made  for  those  cases  where  Search  indi- 
cated p was  less  than  one.  Also  predictions  were  not  made 
for  those  cases  where  p might  have  been  larger  than  one  be- 
cause regression  theory  does  not  explain  how  to  predict  in 
this  case. 

As  seen  in  Tables  II,  III,  and  VI,  the  prediction  in- 
tervals for  the  first  order  autoregressive  model  with  Gen- 
eralized Least  Squares  showed  only  a slight  improvement  over 
the  per  cent  change  regressions.  The  first  order  model  with 

64 


\ 


1 


GOR/SM/76D- 9 


t 


! 

J 


Table  VI 

Generalized  Least  Squares  Prediction 
of  Developmental  Engine  Index 


0 = .95 
GNP^  only 

p = .71 
GNP,  8 

PM? 

1976 

184. 02±  4.60 

193. 19±  4.15 

1977 

208.00+  6.03 

206. 89±  5.26 

1978 

225.66 

223.60 

1979 

245.96 

243.35 

1980 

264.20 

259.73 

1981 

278.50 

273.62 

1982 

297 . 29 

291.89 

1983 

316.54 

310.57 

1984 

329.80 

323.53 

1985 

339 . 4S±18 . 22 

333 . 07±11 . 52 

Prediction  intervals  = 2 op 

where  ap  is  the  standard  deviation  of  the 
prediction 


Generalized  Least  Squares  is,  however,  not  recommended  over 
the  percent  change  technique  for  three  reasons.  As  will  be 
discussed  below,  the  first  reason  is  that  there  is  uncertain-' 
ty  in  the  assumption  of  the  first  autoregressive  model.  The 
second  reason  is  that  the  first  order  autoregressive  model 
with  Generalized  Least  Squares  is  much  more  difficult  to  use. 
The  third  reason  is  that  the  first  order  autoregressive  model 
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would  be  used  for  only  one  of  the  six  indices. 

First  Differences 

The  results  of  the  previous  section  all  indicate  that  p 
is  close  to  one.  When  p is  close  to  one,  first  differences 
is  recommended.  First  differences  was  tried  on  the  six  com- 
posite indices.  However,  the  results  were  judged  inferior 
*,  to  the  per  cent  change  results  for  two  reasons:  per  cent 

changes  had  better  residual  patterns  (see  Chapter  III),  and 

the  first  differences  regression  had  a significant  constant 

I 

term.  The  constant  term  was  statistically  significant  at  the 
95%  confidence  level  for  all  but  the  Development  Engine  equa- 
tion. If  the  assumption  of  a first  order  autoregressive 
model  is  correct  and  if  p is  close  to  one,  then  the  constant 
term  in  the  first  difference  regression  should  not  be  signif- 

t 

icant.  The  significance  of  the  constant  term  indicates  that 
the  assumption  of  the  first  order  regressive  model  with  p 
close  to  one  for  the  model  specified  with  one  independent 
variable,  GNP^,  is  not  valid. 

Economic  principles  reveal  why  first  differences  could 
be  an  inferior  model  compared  to  per  cent  changes.  First 
differences  considers  only  the  change  in  an  index,  not  the 
relative  change.  This  is  fine  when  two  indices  maintain  the 
same  relative  position  to  each  other,  but.  when  one  index  begins 

) 

to  grow  or  decline  much  faster  than  the  other,  a simple 
change  may  no  longer  be  relevant.  It  1958  the  GNP  price  de- 
flator was  66.06  and  the  Development  Airframe  was  64.44,  but 
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in  1975  the  GNP  price  deflator  was  126.37  and  the  Develop- 
ment Airframe  Index  was  176.94.  Relative  to  the  GNP  price 
deflator,  a change  of  say  10  in  the  Development  Airframe  In- 
dex would  be  less  significant  in  1975  than  in  1958.  Per 
cent  changes  would  recognize  this  fact,  first  differences 
would  not. 

The  Primary  Metals  Index,  the  Electrical  Machinery  Index, 
the  national  unemployment  rate,  and  the  lagged  value  of  the 
money  supply  were  all  tried  as  additional  variables  in  the 
per  cent  change  regressions.  None  of  these  variables  were 
recommended  as  additions  to  the  per  cent  change  model  either 
because  they  were  not  statistically  significant  at  the  95% 
confidence  level,  or  that  the  signs  of  their  coefficients 
did  not  make  economic  sense. 

Conclusions  and  Recommendations 

Simply  said,  when  starting  with  the  composite  indices, 
no  technique  was  found  to  be  superior  to  per  cent  change  re- 
gressions. However,  quite  different  results  might  be  obtained 
if  one  started  with  the  subindex  data.  A quick  look  at  the 
subindex  data  shows  three  sets  of  patterns.  One  set  consists 
of  the  labor  indices,  one  set  consists  of  the  material  in- 
dices, and  one  set  consists  of  the  overhead  indices.  It  is 
quite  probable  that  one  could  find  additional  variables  that 
made  economic  sense  and  were  also  statistically  significant 
for  one  of  the  sets  of  indices. 
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VII.  The  Wharton  Predictors 
Errors  in  the  Predictors 

* 

Up  to  this  point  the  predictors  used  in  the  regressions 
equations  were  assumed  to  be  without  error.  The  future,  how- 
ever, is  never  really  known  without  error,  and  the  predic- 
tions used  in  this  thesis  are  no  exception.  This  chapter 
■ will  discuss  the  errors  in  the  predictors  and  how  these  er- 

rors affect  the  predicted  values  of  the  six  aeronautical 
price  indices. 

\ 

The  predictors  used  by  ASD  and  used  in  this  thesis-  ne 
from  the  Wharton  Mark  IV  long  range  econometric  forecast inc 
model.  The  Wharton  model  consists  of  about  1000  interrelated 
equations  which  attempt  to  describe  the  US  economy.  The  mod- 
el is  not  a closed  mathematical  system,  rather  it  needs  exo- 

I 

genous  inputs.  Each  quarter  a group  of  businessmen,  econo- 
mists, and  other  experts  meet  at  Wharton  to  help  determine 
the  subjective  inputs  of  the  mathematical  model.  Besides 
^ business  data,  estimates  of  national  policy  variables  under 

the  influence  of  the  Federal  Reserve,  the  President,  and  the 

I 

Congress,  are  placed  into  the  model.  Wharton  does  allow  for 
multiple  solutions  depending  upon  the  user's  preference.  For 
example,  Wharton  provides  estimates  of  the  economy  for  a 
tight  money  or  easy  money  scenario.  The  fact  remains, 
though,  that  the  Wharton  estimates  depend  on  both  mathemati- 
cal equations  and  subjective  judgements  (Ref  26:  Ref  27). 

There  are  models  that  do  not  use  subjective  inputs. 
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Brush,  for  instance,  proposes  a prediction  for  inflation  that 
uses  no  (or  at  least  very  little)  subjective  inputs  (Ref  28). 
Brush's  model  may  give  good  results  for  two  years  into  the 
future.  This  is  because  all  economic  policy  decisions,  es- 
pecially those  of  Congress,  have  a considerable  lag  time  be- 
fore implementation.  Furthermore,  implemented  policy  takes 
time  before  its  effect  is  felt.  Unfortunately,  ASD  has  a 
need  for  estimates  up  to  ten  years  in  the  future.  The  need 
for  long  term  predictions  and  a desire  for  consistency  in 
reporting,  caused  ASD  to  choose  Wharton's  predictors  over 
Brush ' s . 

How  good  are  the  Wharton  forecasts?  That  is  hard  to 
say.  Only  one  study  could  be  located  that  objectively  eval- 
uated Wharton's  predictors  (Ref  29).  It -compared  a number 
of  economic  forecasting  models.  Wharton's  Mark  IV  was  not 
included,  but  rather  the  older  Mark  III  quarterly  version 
was  tested.  However,  the  Wharton  model  did  show  the  lowest 
error  of  all  the  models  tested  over  a two  year  period.  The 
errors  for  the  Wharton  GNP  price  deflator  for  two  years  in 
the  future  were  less  than  3 percentage  points  of  the  1958 
value . 

How  the  Errors  Affect  the  Price  Index  Predictions 

How  do  the  errors  in  the  Wharton  predictors  affect  the 
errors  in  the  aeronautical  price  indices?  Figs.  12  and  13 
will  help  explain  the  effect.  In  this  example  the  Develop- 
ment Engine  index  will  be  used.  The  relationship  between 
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the  Development  Engine  index  (Dev  Engine)  and  the  GNP  price 
deflator  (GNP^)  is  represented  by  line  1 in  Fig.  11.  For 
any  GNP^,  e.g.  point  A,  a corresponding  estimate  of  Dev  Eg- 
gine,  point  B,  can  be  determined.  But  what  if  point  A is 
not  known  exactly  but  rather,  because  of  error,  is  known  to 
be  in  the  interval  CD.  Then  if  there  were  no  error  in  the 
Dev  Engine  relationship,  the  predictors  of  Dev  Engine  would 
be  in  interval  FE.  If  FE  is  larger  than  CD  (assuming  equal 
scales)  then  this  is  the  penalty  one  pays  when  he  wants  to 
predict  a specific  index,  not  just  GNP^ . 

The  relationship  between  DAF  and  GNP^  is  not  known  with- 
out error.  Prediction  intervals  for  their  relationship  is 
represented  by  lines  2 and  3 in  Fig.  12.  Now,  the  interval 
CD  results  in  interval  HG,  which  is  better  than  FE.  The 
phenomena  represented  here  in  Figs.  11  and  12,  is  not  quite 
mathematically  exact,  but  it  does  represent  what  really  hap- 
pens . 

Adding  sample  numbers  will  illustrate  the  magnitude  of 
the  errors.  From  Table  III,  the  predicted  error  in  1977  for 
the  Development  Engine  Index  is  6.67.  The  root  mean  square 
error  of  the  Wharton  predictor  is  approximately  3%  of  the 
1958  value  of  66.06,  which  equals  1.98.  If  one  considers  the 
root  mean  square  error  as  an  approximation  to  the  standard 
deviation  of  the  error,  then  two  times  the  root  mean  square 
error  would  yield  95%  confidence  bounds,  if  the  error  was 
normally  distributed.  Thus  the  error  in  GNP^  would  be  3.96. 
The  regression  line  for  the  indices  has  an  equivalent  slope 


GOR/SM/76D- 9 


of  close  to  1,  so  the  total  error  for  the  Development  Engine 
Index  in  1977  would  be  ±1 . 98  + 6 . 67  = ±8 . 65  . This  is  certainly 
only  an  estimate,  and  is  not  mathematically  exact,  but  it 
does  give  an  idea  of  the  magnitude  of  the  possible  predic- 
tion errors. 
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VIII.  Summary,  Findings , and  Recommended 
Areas  for  Future  Study 

Summary 

The  thesis  began  with  a description  of  the  annual  cost 
escalation  reports  prepared  by  the  Aeronautical  Systems  Di- 
vision (ASD)  of  the  United  States  Air  Force  Systems  Command. 
The  reports  contain  predictions  for  selected  aircraft  price 
indices  for  up  to  ten  years  in  the  future.  The  thesis  then 
recalculated  much  of  the  historical  price  index  data  used  by 
ASD  in  their  reports.  Next  the  price  indices  were  repredicted 
using  the  same  methodology  as  used  by  ASD  in  their  latest  re- 
port. Prediction  intervals  were  also  included.  The  thesis 
then  examined  the  question  of  whether  it  was  better  to  aggre- 
gate the  historical  data  first  and  then  predict,  or  predict 
subsets  of  the  historical  data  and  then  aggregate  the  predic- 
tions. Finally,  advanced  statistical  techniques  were  tried 
in  an  attempt  to  improve  the  predictions. 

Findings 

The  significant  findings  of  this  thesis  can  be  separated 
into  four  different  areas:  regression  theory,  historical 
data,  regression  techniques,  and  prediction  intervals. 

Regression  Theory 

The  use  of  the  Durbin  Watson  statistic,  D,  is  incorrect- 
ly explained  in  some  textbooks.  The  significance  levels 
given  in  the  Durbin  Watson  tables  only  apply  to  cases  where 
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D > D^ , i.e.,  when  the  hypothesis  of  no  autocorrelation  is 
rejected . 

The  Cochrane-Orcutt  iterative  method  for  estimating  the 
parameters  of  a first  order  autoregressive  model  is  frequently 
illustrated  without  the  first  row  of  its  transformation. 

The  use  of  Cochrane-Orcutt  without  the  first  row,  in  cases 
where  autocorrelation  has  been  going  on  before  the  first 
data  point,  reduces  the  efficiency  of  the  technique. 

The  Search  technique  developed  in  Chapters  III  and  VI, 
yields  maximum  likelihood  estimates  of  the  parameters  of  the 
first  order  autoregressive  model.  The  Search  technique  is 
recommended  by  the  author  as  the  most  accurate  technique  for 
estimating  the  parameters  of  the  first  order  autoregressive 
model . 

Historical  Data 

Almost  all  of  the  historical  data  used  by  ASD  in  their 
latest  report  was  examined  and  recalculated.  There  were  only 
a few  minor  differences  between  the  recalculated  data  and  the 
data  listed  in  the  ASD  HOC  report.  The  most  significant 
difference  was  in  the  calculation  of  the  Airframe  Raw  Material 
subindex.  The  difference  was  apparently  due  to  a difference 
in  the  ASD  calculations  when  the  Bureau  of  Labor  Statistics 
switched  from  a 1957-1959  base  to  a 1967  base.  ASD  has  re- 
calculated the  Airframe  Raw  Material  Subindex  values  for  use 
in  future  reports. 
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Regression  Techniques 

Given  the  independent  variables  listed  in  the  latest 
ASD  HOC  report,  the  GNP  Price  Deflator,  the  Primary  Metals 
index,  and  the  Electrical  Machinery  index,  the  following 
conclusion  was  made.  The  aggregate  per  cent  change  method 
is  preferred  slightly  to  the  per  cent  change  subindex  method. 

When  starting  with  the  composite  data,  the  per  cent 
change  regression  technique  was  found  to  be  superior  to  re- 
gression with  nonlinear  variables,  multivariate  regression, 
Generalized  Least  Squares,  and  first  differences. 

Prediction  Intervals 

Using  aggregate  per  cent  change  regression,  the  predic- 
tion intervals  for  the  six  composite  aircraft  price  indices, 
given  no  error  in  the  Wharton  predictions,  varied  widely. 
However,  they  averaged  approximately  1.81  for  the  first 
year,  2.51  for  the  second  year,  and  4.71  for  the  tenth  year 
of  the  predicted  values.  Errors  in  the  Wharton  predictions 
would  be  expected  to  raise  the  error  in  the  predictions  to 
a total  of  3.5%  for  the  second  year. 

Recommended  Areas  for  Future  Study 

The  historical  data  is  an  important  ea  for  future 
study.  The  most  advanced  regression  techniques  and  the  most 
accurate  predictions  of  the  independent  variables  will  be 
useless  if  the  historical  data  does  not  accurately  represent 
the  true  historical  price  changes.  The  historical  price  in- 
dices were  recalculated  in  this  thesis,  but  the  ASD  metho- 
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dology  used  to  calculate  the  indices  was  assumed  correct. 

The  methodology  should  be  carefully  examined. 

The  composite  indices  were  used  as  a starting  point  when 
investigating  other  possible  regression  techniques.  The  same 
investigation  could  be  done  starting  with  the  subindex  data. 

Finally,  the  variance  of  the  predictions  was  calculated, 
but  the  probability  distribution  of  the  predictions  was  un- 
known. The  normal  distribution  was  assumed  in  order  to  cal- 
culate prediction  intervals.  Monte  Carlo  techniques  could 
be  used  to  construct  sample  probability  distributions  leading 
to  perhaps  more  accurate  prediction  intervals. 
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Appendix  C 

Calculation  of  Prediction  Intervals 
for  Per  Cent  Change  Regressions 

This  appendix  shows  the  theory  used  to  estimate  the  pre- 
diction intervals  that  are  contained  in  Chapter  V.  The  pre- 
dictions are  for  the  aeronautical  indices  that  were  derived 
by  two  methods:  the  subindex  per  cent  change  method,  and  the 
aggregate  per  cent  change  method.  Calculation  of  prediction 
intervals  for  the  subindex  method  required  six  steps,  while 
calculation  of  prediction  intervals  for  the  aggregate  method 
required  only  four  steps. 

Prediction  Intervals  for  the  Subindex  Method 


Step  1^  Check  Assumptions  of  Ordinary  Least  Squares 

The  assumptions  of  Ordinary  Least  Squares,  discussed  in 
Chapter  III,  were  checked  for  each  regression  equation.  This 


ensured  that  the  errors  within  each  regression  were  indepen- 
dent and  were  normally  distributed. 


i 

■ 


Step  2 Calculate  Variances  for  the  Per  Cent  Change  Prediction 
The  Wharton  perdictors  were  substituted  in  to  the  re- 
gression equations  to  yield  predictions  for  the  per  cent 


change  of  the  subindices.  The  variances  for  these  per  cent 


i 


change  predictions  were  easily  calculated  using  the  following 
formula : 


KL 
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where  Var(Aly)  is  the  variance  of  the  per  cent  change  for 
2 2 

year  y.  a is  estimated  by  S , and  is  the  vector  of 

predictors  for  year  y (Ref  5:135).  The  results  for  the  en- 
gineering labor  subindex  were: 


Engineering  Labor  = 1.582  + .8286  A%GNP( 
A% . 


'76 
Var(A%76) 

A I 


= 6.4937 


'll 
Var (A%  ??) 


4.25339 

6.6707 

4.287 


(65) 

(66) 

(67) 

(68) 
(69) 


Step  3_  Calculate  Variances  of  the  Subindex  Predictions 

The  predictions  are  easily  calculated  using  Eqs  (50) 
and  (51).  The  variances  are  more  difficult  to  calculate. 
The  following  is  an  example  of  the  calculations  for  the  en- 
gineering labor  subindex  along  with  a brief  explanation: 


A% 


1976  =(W-+  (1975)  = (.6^~  + 1)  (171.81)  = 182.97 


100 


(70) 


-1975, 


.171.81. 


Var  (1976)  = (^f)  Var(A%y6)  = ( l00~)  (4.25339)  = 12.558 


100 
A%. 


(71) 


1977  = (T^07;!  + ^ (1976)  = (6  + 1)  (182.97)  = 195.17 

(72) 


A% 


Var(1977)  =Var[(-j-^  + 1)  (1976)  ]=  y^-Var  [ (A*  ?y)  (1976)  ] + 


Var  1976  = [ (1976)  2 Var(AI?7)  + (A%?6)2  Var  (1976)  + 


Var  (1976)  Var(A%?6)]  + Var(1976)  = [ (182 . 97)  2 (4 . 287)  + 


(6.6707)^(12.555)  + 112 . 555) (4 . 287) 1 + 12.55  = 26.968 

U 3) 
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where  1976  means  the  1976  prediction  of  the  subindex.  Eq 
(73)  uses  the  assumption  that  the  1975  index  is  known  with- 
out error  and  that  the  predictors  are  independent  of  each 
other.  Eq  (73)  is  based  on  the  following  relationship: 

if  X and  Y are  independent  random  variables  then 
Var(XY)  = E(X)2  Var(X)  + E (Y) 2 Var(Y)  + Var(X)  Var(Y)  (74) 

where  E(X)  is  the  expected  value  of  X. 

Step  £ Estimate  the  Covariance  Between  the  Subindices 

The  covariances  are  needed  for  Step  5.  They  are  esti- 
mated using  the  sample  correlation  matrix  of  the  residuals 
of  the  subindex  regressions  adjusted  each  year  for  the  vari- 
ances of  the  subindex  predictions.  The  following  is  an  ex- 
ample: 

sample  correlation  between  subindex  5 (Avionics  Devel- 
opment Labor)  and  subindex  9 (Avionics  Raw  Material) 


« j359  = -.  2444  (75) 
Variance  for  subindex  5,  1976  = ct^2  = 4.027  (76) 
Variance  for  subindex  9,  1976  = a^2  = 7.809  (77) 
Covariance  between  5 and  9 = °59  = 

-.  2444  (/47 027)  (/7."809)  - -1.378  (78) 


Step  _5  Determine  Variances  of  Composite  Indices 

The  composite  indices  are  simply  a linear  combination 
of  subindices  as  shown  in  Eqs  (19)  through  (26).  The  vari- 
ances of  a linear  combination  of  interdependent  random  vari- 
ables is  calculated  as  follows: 
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if  Y = aX  + bX 

then  Var(Y)  = a2a2x  + 2aboxz  = b2a2z  (79) 

where  X,  Y,  Z are  random  variables,  a,  b are  constants  and 
2 

ox  is  the  variance  of  X,  and  oxz  is  the  covariance  between 
X and  Z . 

So  for  the  Avionics  Development  Index 

Var  (1976)  = (. 9693)  2 (4 . 027)  + (.1)  2(7.889)  + 

2(.l)  (.9693)  (-1.378)  = 3.595  (80) 

Step  6 Calculate  Prediction  Intervals 

Step  1 started  with  normally  distributed  random  errors. 
The  predictions  in  Step  2 were  still  normally  distributed, 
but  in  Step  3 the  predictions  were  no  longer  normally  dis- 
tributed. In  Step  4 only  estimates  of  the  covariances  were 
available.  So  the  variances  in  Step  5 are  only  estimates 
and  they  are  for  unknown  probability  distributions.  Still 
it  is  felt  that  two  times  the  standard  deviation  would  be 
close  to  a 95%  prediction  interval  (Ref  30:8).  Thus  for  the 
Avionics  Development  Index  a 95%  prediction  internal  for 
1976  is 

= 180.27±2/3.595 

= 180 . 27±  3 . 7 (81) 

Prediction  Intervals  for  the  Aggregate  Method 

Step  1^  Check  Assumptions  of  Ordinary  Least  Squares 

This  is  exactly  the  same  as  for  the  subindex  method, 
except  now  the  six  aggregate  regressions  are  Checked. 


G0R/SM/76D- 9 

Step  2^  Calculate  Variances  for  Index  Predictions 

As  in  Step  2 of  the  subindex  method,  the  Wharton  pre- 
dictors were  substituted  in  the  regression  equations  to  yield 
the  prediction  for  the  indices.  The  variances  were  again 
easily  calculated.  For  example,  for  the  Development  Air- 
frame Index: 


[ * 

N 


i 


i 

* 


i 

' 


I * 


; » 


j 


4 


A%  Development  Airframe 
A%76  = 7.2455 
Var(A%?6)  = .42362 
A %77  = 7.3640 
Var  (A%  77)  = .42701 


(DAF)  = 3.959  + .5544  A%  GNPd 

(82) 


Step  Calculate  Variances  of  Index  Predictions 

The  same  mathematical  theory  that  was  used  in  Step  3 of 
the  subindex  method  is  used  in  this  step.  Examples  for  the 
Development  Airframe  Index  are: 


1976  = 176.9454  (^^  + 1)  = 189.765  (83) 

2 

Var  (1976)  = C176^454)  (.42362)  = 1.326  (84) 

1977  = 189.765  (7j^4  + 1)  = 203.739  (85) 

Var  (1977*)  = [ (189 . 765)  2 (.  42701)  + (7 . 364)  2 (1 . 326)  + 

(1.326)  (.42701)]  + 1.326  = 2.87  (86) 


Step  4 Calculate  Prediction  Intervals 

Since  the  data  was  already  aggregated,  covariances  are 
not  needed.  Through  Step  2 the  variances  were  normally  dis- 


91 


tributed.  However,  the  variances  from  Step  3 are  not  normal- 
ly distributed,  but  it  is  still  felt  that  two  times  the  stan- 
dard deviation  will  yield  good  estimates  for  951  prediction 
intervals.  For  example,  the  prediction  interval  for  the  1976 
Development  Airframe  index  is: 


189.765  ± 2/1.326 

189.765  ± 2.3  (87) 


GOR/SM/76D-9 


Appendix  D 


Flow  Charts  and  Computer  Codes 


GOR/SM/76D- 9 


Durbin  Two  Step  Computer  Code  (Ref  6) 


Read  in  subindices 

own l TAR 


DIMENSION  98X4? 

SET  ENGINEERING  LABOR  IN  1 

77.52  79.0?  82.55  84.23  83.72  86.41  91.44  94.46  95.64  100, 
107.21  115.1  125.34  130.03  137.08  144.96  157.05  171.81 
SET  MANIIF ACT  LABOR  IN  COL  2 
71.92  75.64  77.65  79.66  82. ?3  84.53  85.96 

104.30  111.75  119.48  126.65  138.97  146.99 
SET  AC  PARTS  AND  EQUIP  IN  COL  3 
72.84  76.1?  78.80  80.60  83.58  86.27  88.96 
105.37  112.24  119.10  124.48  132.54  139.70 
SFT  AC  ENGINE  LABOR  IN  COI  4 
73.39  77.19  79.82  82.16  85.09  87.43  90.35 
106.72  113.16  119.88  128.07  139.18  147.95 
SET  AVIONIC  DEV  LABOR  IN  COL  5 
68.82  71.96  74.71  78.04  81.18  84.51  87.65 
106.86  113.33  120.00  126.67  132.94  140.00 
SET  AVIONIC  PROD  LABOR  IN  COL  6 
71  .13  74  . ?r  76.64  80.05  82.41  85.56  88.19 

106.30  112.60  119,42  125.46  131.76  140.42 
SFT  AT  RAW  MAT  IN  COL  7 

100.31  96,58  99.67  99.77  98.96  95.49  96.61  96.90  97.89  100. 
101.70  107.92  114.28  113.60  113.32  116,32  148.53  165.92 
SET  ENGINE  RAW  IN  COL  8 

112.62  108.44  106.58  107.10  105.80  100.89  99.02  99.54  98.31 
100.  102.88  107.51  118.82  122.17  123.36  135.27  165.38  200.18 
SFT  AVIONIC  RAW  MAT  IN  COl  9 

105.4  104.1  102.8  100.9  100.6  100.  97,0  96.  98.4  100. 

101.3  103.4  105.2  108.2  109.6  110.3  118.1  127.8 
SET  AVIONIC  DEV  LABOR  IN  COL  5 

68.82  71.96  74.71  78.04  81.18  84.51  87,65  90.39  93.14  100. 
106.86  113.33  120.00  126.67  132.94  140.00  148.69  162.12 
SET  AVIONIC  PROD  IAROP  IN  COL  6 

71.13  74.28  76.64  80.05  82.41  85.56  88.19  91.08  93.70  100. 

106.30  112.60  119.42  125.46  131.76  140.42  148,92  163.04 
SET  AF  PAW  MAT  IN  COL  7 

100.31  96.58  99.67  99.77  98.96  95.49  96.61  96.90  97.89  100. 
101.70  107.92  114.28  113.60  113.32  116,32  148.53  165.92 
SET  ENGINE  RAW  IN  COL  8 

112.62  108.44  \06.58  107,10  105.80  100.89  99.02  99.54  98.31 
100.  102.88  107.51  118.82  122.17  123.36  135.27  165.38  200.18 
SET  AVIONIC  RAW  MAT  IN  COl  9 

105.4  104,1  102.8  100.9  10U.6  100.  97.0  96.  98.4  100. 

101,3  103.4  105.2  108.2  109.6  110.3  118.1  127.8 


90.26  95.70  100. 
159.60  177.65 

91.94  95.8?  100. 
149.55  164.78 

92.69  97.08  100. 
158.77  176.32 

90.39  93.14  100. 

148.69  162.12 

91.08  93.70  100. 
148.92  163.04 


95 
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SFT  SFVFN  PT  TuREF  0 VF  RHEA  D IN  COL  10 

*53.0-1  56  . PI  61.07  65.52  70  . .11  75.44  80.95  86.86  91.20  10  0 . 
107.10  115.13  123.54  132.86  142.21  152.61  163.75  175.71 
SET  STVEN  PT  FIGHT  OVERHEAD  IN  11 

50.87  54.81  5P.11  63.72  68.69  74.05  79.83  86.05  92.76  100. 
107.8  116.21  125.27  135.04  145.58  156.93  169.17  182.37 
SFT  AVIONICS  PROD  OH  IN  11 

63.17  66.84  60.88  73.96  77.15  81.16  84.77  88.71  92.47  100. 
107.71  115.61  124.24  132.26  140.75  151,98  163.32  181.18 
SFT  AVIONICS  DFV  OH  IN  12 

61.12  64.76  68.12  72.1  76.  80.17  84.25  88.04  91.92  100. 
108.28  116.36  124.85  133.54  142.  151.52  163.07  180.16 

Aggregate  data  and  Regress 


i 


i 

I 

i 

i 


MIILT  7 BY  .4  PUT  IN  14 

MULT  3 BY  .225  BUI  T BY  1.  ADD  14  PUT  IN  14 

MULT  10  BY  .375  1.  1.4  14 

MULT  2 RY  .1428  PUT  IN  19 

MULT  1 BY  .2125  1.  19  19 

MULT  14  RY  .1924  1.  19  19 

MIILT  11  BY  .4526  1 . 19  19 

MULT  2 BY  .2010  PUT  IN  20 

MULT  1 RY  .0650  1.  20  20 

MULT  14  BY  ,3040  1.2020 

MULT  11  BY  .4300  1.  20  20 

MULT  8 BY  .4  PUT  IN  15 

MULT  4 BY  .6  BY  l.  15  15 

MULT  4 BY  .2701  PUT  IN  21 

MULT  1 BY  .3415  1.  21  21 

MULT  8 BY  .1942  1.  21  21 

MULT  15  RY  .1942  1.  21  21 

MULT  4 BY  .4184  PUT  IN  22 

MULT  1 BY  .0416  1.  22  22 

MULT  8 BY  .270(1  1.  22  22 

MULT  15  BY  .2700  1.  22  22 

MIILT  5 BY  .15  PUT  IN  21 

MULT  9 BY  .1  1.  23  23 

MULT  12  BY  .55  1.  23  23 

MULT  6 BY  .3  PUT  IN  24 

MULT  9 BY  .2  1.  24  24 

MULT  13  BY  .5  1 . 24  24 

MOVE  1.19  18X6  TO  1,1 

MOVE  2.1  17X9  TO  1,10 

SFT  GNP  DFFlATnR  IN  COL  30 

66.06  67.52  68.67  69.28  70.55  71.59  72.71  74.32  76.76  79.02 
82.57  86.72  91.36  96.02  100.  105.92  116.2  126.37 
DFFINE  1.0  IN  COL  28 
ERASE  7***27 
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MOVE  2,  1 17X6  TO  1.7 
MOVE  2,3(1  17X1  TO  1,31 
SPACE  3 

MPR 1 NT  1,1  18X12 

MPRJNT  1,30  18X2 
RESET  17 
REG  1 N 

iriT  Y IN  10  WTS  1.0  4 VAR  IN  2 8,31,30,4  C 33 
2SPACE  3 

3MPR1NT  1,33  18X1 

4 I NCRFMENT  INSTRUCT  1 BY  1 , 0 . , 0 , 0 , 0 , 0 , 1 , 0 
FINISH 

PERFORM  1 TO  4 6 T/Mf-S 

STOP 
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Search  Computer  Code  (Ref  6) 
Read  in  subindices 
Omni TAB 


DIMENSION  28X4? 

SFT  ENGINFFRINC  LAGOP  IN  1 

77.52  79.02  82.55  84.23  83.72  86.41  91.44  94.46  95.64  100. 
107.21  115.1  125.34  130.03  137.08  144.96  157.05  171.81 
SFT  MANUFACT  LABOR  IN  COL  2 

71.92  75.64  77.65  79.66  82.23  84.53  85.96  90.26  95.70  100. 

104.30  111.75  119.48  126.65  138.97  146.99  159.60  177.65 
SFT  AC  PARTS  AND  FQIJIP  IN  COL  3 

72.84  76.12  78.80  80.60  83.58  86.27  88.96  91.94  95.82  100. 
105.37  112.24  119.10  124.48  132.54  139.70  149.55  164.78 
SFT  AC  ENGINE  LABOR  IN  COI  4 

73.39  77.19  79.82  82.16  85.09  87.43  90.35  92.69  97.08  100. 
106.72  113.16  119.88  128.07  139.18  147.95  158.77  176.32 
SET  AVIONIC  DEV  LABOR  IN  COL  5 

68.82  71.9^  74.71  78.04  81.18  84.51  87.65  90.39  93.14  100. 
106.86  113.33  120.00  126.67  132.94  140.00  148.69  162.12 
SET  AVIONIC  PROD  LABOR  IN  COL  6 

71.13  74.28  76.64  80.05  82.41  85.56  88.19  91.08  93.70  100, 

106.30  112.60  119.42  125.46  131.76  140.42  148.92  163.04 
SET  AF  RAW  MAT  IN  COL  7 

100.31  96.58  99.67  99.77  98.96  95.49  96,61  96.90  97.89  100. 
101.70  107.92  114.28  113.60  113.32  116.32  148,53  165.92 
SFT  ENGINF  RAW  IN  COL  8 

112,62  108.44  106.58  107.10  105.80  100.89  99.02  99.54  98.31 
100.  102.88  107.51  118.82  122.17  123.36  135.27  165.38  200.18 
SFT  AVIONIC  RAW  MAT  IN  COL  9 

105.4  104.1  102.8  100.9  100.6  100.  97.0  96.  98.4  100, 

101.3  103.4  105.2  10R.2  109,6  110.3  118.1  127.8 
SFT  PRIMARY  METAL  DEFLATOR  IN  COL  41 

87.26  88.57  91.54  90.05  92.24  91.71  92.85  93.46  95.20  100. 
105.85  109.6  113.88  120.34  127.58  126.61  144.34  167.97 
178.59  192.03  ?06.36  220,72  235.74  248.38  264.89  281.58 
293.47  302.53 


SET  ELECT  mFTAI  DEFLATOR  IN  COL  42 

113.25  114.27  112.  111.1  107.25  104.53  102.15  99.21  97.73  100. 

102.72  98.64  103.17  105.21  105.1  106.23  113.25  124.58 
135.90  147.23  158.55  169.88  181.20  192.53  203.85 
215.18  226.5  237.83 

SFT  SFVFN  PT  THREE  OVERHEAD  IN  COL  10 

53.04  56.91  61.07  65.52  70.31  75,44  80.95  86.86  93.20  100. 
107,30  115.13  123.54  132.56  142.23  152.61  163.75  175.71 
SET  SEVEN  PT  FIGHT  OVERHEAD  IN  11 

50.87  54.83  59.11  63.7?  68.69  74.05  79.83  86.05  92.76  100. 
107.8  116,21  125.27  135.04  145.58  156,93  169.17  182,37 
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SFT  AVIONICS  PROP  OH  IN  13 

63.17  66. #4  69.88  73.06  77.15  81.16  84.77  88.71  9?. 47  100. 
107.71  115.61  124.24  132. ?6  140.75  1 51.98  ]63.3?  181.18 
SFT  AVIONICS  OP  V OH  IN  12 

61.12  64.76  68.12  72.10  76.  80.17  84.25  80.04  91.92  100. 
108.28  116.36  124.85  133.54  142.  151.52  163.07  180.16 


Aggregate  data 


MULT  7 BY  .4  PUT  IN  14 

MULT  3 BY  .225  BULT  BY  1.  ADO  14  PUT  IN  14 

MULT  10  BY  .375  1.  14  14 

MULT  2 BY  .1425  PUT  IN  19 

MULT  1 BY  .2125  1.  19  19 

MULT  14  BY  .1924  1.  19  19 

MULT  11  BY  .4526  1.  19  19 

MULT  2 BY  .2010  PUT  IN  20 

MULT  1 BY  .0650  1.  20  20 

MULT  14  BY  .3040  1.  20  20 

MULT  11  BY  .4300  1 . 20  20 

MULT  8 BY  .4  PUT  IN  15 

MULT  4 BY  .6  BY  1.  15  15 

MULT  4 BY  .2701  PUT  IN  21 

MULT  1 BY  .3415  1.  21  21 

MULT  8 BY  .1942  1.  21  21 

MULT  15  BY  .1942  1.  21  21 

MULT  4 BY  .4184  PUT  IN  22 

MULT  1 BY  .0416  1.  22  22 

MULT  8 BY  .2700  1.  22  22 

MULT  15  BY  .2700  1.  22  22 

MULT  5 BY  .35  PUT  IN  23 

MULT  9 BY  .1  1.  23  23 

MULT  12  BY  .55  1.  23  23 

MULT  6 BY  .3  PUT  IN  24 

MULT  9 BY  .2  1.  24  24 

MULT  13  BY  .5  1.  24  24 

MOVE  1,19  18X6  TO  1,1 

ERASF.  7*  * # 24 

SET  UNEMPLOYMENT  IN  COL  18 


6.8  5.5  5.5  6.7  5.5  5.7 
3.6  3.5  4.9  5.9  5.6  4.9 
SET  LAGOEB  M ONE  IN  COL 
135.9  141.1  143.4  144.2 
156.5  163.7  171.3  175.4 
221.4  235.3  255.8  271.5 


5.2  4.5  3.8  3.8 
5.6  8.5 
17 

148.7  150.9 
186.9  201,7  208.7 
284.4 
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RESET  28 

define  t.o  in  ool2b 

SFT  ONP  DFFLATOR  IN  COL  30 

66,06  67.5?  6 8.67  69. ?8  70.55  71.50  72.71  70.3?  76.76  79.02 
8?.57  86.7?  91.36  96.0?  100.  105.92  116.2  126.37 

133.86  142.08  152.48  164.45  175.20  183.62  194.70  206.06 

213.87  219.55 
SPACE  3 

MOVE  1.5  18X1  TO  1,1  (INDEX  FIVE) 

! " NOTEUWORKSHFET 

MPRINT  1,1  18X6 
MPRINT  1,17  18X2 
MPRINT  1,30  18X1 

NOTE  1 % RHO,  COFFE , DURBIN  WATSON 
DEFINE  1.0  IN  COL  13 
: RFSET  18 

BEGIN 


Initialize  Rho , p 


1ADEFINF  1,8  18X1  TO  B E 0.80  ( RHO) 
2 SOU A RE  8 PUT  IN  9 
3 SUB  TRACT  9 FM  28  PUT  IN  10 
4S0RT  OF  10  IN  U 


' Transform  composite  Price  Index 

1 

5AMULT  1,1  18X1  BY  11  PUT  IN  1,7 
\ 6AMIJLT  1,1  17X1  BY  8 PUT  IN  2,7 

7SUBTRACT  7 FM  1 PUT  IN  20 


Transform  GNP 


: * 8AMIILT  1,30  18X1  BY  11  PUT  IN  1.31 

9AMULT  1,30  17X1  BY  8 PUT  IN  2,31 
j - 10SUBTRACT  31  FM  30  PUT  IN  ?1 

[ i 
\ 

Transform  constant 

% 

11AMULT  1,13  18X1  BY  11  PUT  IN  1.12 
12AMULT  1,13  17X1  BY  8 PUT  IN  2,1? 
13SUB TRACT  1?  FM  13  PUT  IN  22 
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Regress 


17SFIT  Y IN  20  WTS  1.0  2 VR  !N  ??.,  21  C IN  3J  « IN  34 


Calculate  Durbin  Watson  Statistic 


18M0VE  2,34  17X1  TO  1,35 

19SUBTRACT  34  FM  35  PUT  IN  35 

20SOUARE  34  PUT  IN  36 

21SQUARF  35  PUT  IN  37 

22SUM  36  ROW  1 TO  18  PUT  IN  38 

23SUM  37  ROW  1 TO  17  PUT  IN  39 

24DIVIDE  *1,39*  RY  *1,38*  PUT  IN  40 


Print.  Rho , Coefficients  including  S^,  Durbin  Watson  Statistic 


25SPACE  3 
R6PRINT  NOTF 
?7SPACE  3 
28MPRINT  1,8  1X1 
29SPACE  3 

30MPR1NT  1,33  10X1 


31SPACE  3 

32MPRINT  1,40  1X1  <DW) 


Increment  Rho 


33 1 NCREMENT  INSTRUCT  1 RY  0,0,0,0,0.01 
FINISH 

PFRFORM  1,33  19  TIMES 


Next  Composite  Index 


RESTORE  1 TO  1,8,18  ,1 ,0.80 
MOVE  1,6  18X1  TO  1,1 
PERFORM  1,33  19  TIMES 
STOP 
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